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Summary. 

We report several results of orbit determinations of comet 
Halley. Some problems which appear are considered, especially in 
regard to the nongravitational forces and the differences between 
the light center and the nucleus of the comet. Improved orbital 
elements have been computed for different assumptions about these 
and other bases. Further results are: a radial light offset 
can most simply to be eliminated by using only the position angles 
to sun of the observations (instead a and 5); due to the decay of 
the comet the nongravitational forces increase by about 1% per 
revolution, the nongravitational forces decrease much slower at large 
heliocentric distances than according to the nongravitational models, 
and observations back to at least 1759 should be used for orbit 
determinations until the comet is on larger distances. The perihelion 
time in A.D. 837 has been determined very precisely and later can be 
used as a check of improved models of the nongravitational forces. 
A backward integration of the comet back to 2300 B.C. has been carried 
out and compared with the observed perihelion times. 



Zusanmenfassung. 

Es wird fiber Ergebnisse der Bahnberechnungen des Halley'schen Kometen 
berichtet. Die auftretenden Schwierigkeiten werden erortert, insbesondere 
in Hinblick auf die nichtgravitativen Krafte und die Differenz zwischen 
dem Lichtschwerpunkt und dem Kern des Kometen. Zu verschiedanen Aimahmen 
fiber diese und andere Grundlagen wurden verbesserte Bahnelemente abgeleitet. 
Weitere Befunde sind: einen radialen light offset kairn man am einfachsten 
eliminieren, indem man nur die Positionswinkel zur Sonne (statt a und 5 ) von 
den Beobachtungen verwendet, durch Aliening des Kometen vergroSern sich die 
nichtgravitativen Parameter um etwa 1% pro Umlauf, die nichtgravitativen 
KrMfte nehmen mit grofler Entfernung erheblich langsamer ab als gemaQ den 
vorhandenen Modellen, und bei Bahnrechnungen sollten deshalb Beobachtungen 
zuriick bis mindestens 1759 verwendet werden, solange der Komet in groBeren 
Distanzen ist. Die Perihelzeit im Jahre 837 wurde sehr genau bestimmt und 
kann zur tJberpnifung verbesserter Modelle fur die nichtgravitativen Krafte 
verwendet werden. Eine Riickrechnung des Kometen bis 2300 v.Chr. wurde 
durchgefiihrt und mit den beobachteten Periheldurchgangen verglichen. 



" : 0*1 town,*. u 

c^^M. . y^H(>u<a&T*ioH- CUIL6JU, \t*tt -6 c^Ufr6&**i'<@'%t-tUfru<ir 



0. Introduction 

After the rediscovery of comet Hal ley in October 1982, the author 
started computations with the aim to derive improved orbital elements 
for the comet (w.Landgraf, 1983a). The present report reviews the 
proper results of these computations, together with comments on the 
bases of them and their uncertainties and difficulties, especially in 
regard to the nongravitational forces. In addition, some other 
aspects are pointed out which should be taken into account in similar 
computations. 

1. The observations 

Reliable Reports refering to comet Halley extend back to 466 B.C. 
The results presented here refer only to the 1607 to 1984 appearences. 
However, they exhibit significant differences to the earliest reports, 
and thus the latter are principially of value for more extended 
investigations of the long-term motion of the comet. 

From 1909 to 1911 about five thousand observations have been reported 
and about half of their most precise ones have been reduced to the 
system of GC and collected to 33 normal places (p.zadunaisky, 1966} . 
These normals has been choosen for our computations and reduced from GC 
to the FKit system. Although the time of brighter magnitude of the 
comet is covered by a larger number of observations, the greater 
observational uncertainty there and consideration of the residuals of the 
normals favours weighting all normal places with one unit weight, which 
a posteriori corresponds to the mean error of 1"2 . In addition, four 
observations recently measured (B.Bowell 1982) were used with unit weight 
in a and a mean weight of 0.6 in 6. 



About 400 observations are published on the apparition of 1835. 
These are, however, of very different quality. The observations 
by j.Maclear (1837) at Cape, for example, immediately show a 
scattering of some I s in subsequent right ascensions, while in those 
by j.F.w.Herrschel (1837) the observation times are given only to one 
minute. Other observations by different observers are either from 
rather poor observation series, or have shown both before and after 
re-reduction large residuals. Some of the most favourable observation 
series are these by j.Encke (1838) at Berlin and j.Lamont (1843) at 
Bogenhausen. From the first observer, after removing each three right 
asecnsions and declinations with much large residuals of the 25 
observations which were re-reduceable using modern comparison star 
positions, there remains 44 measurements between 1835 September 18 
and 1836 March 19 with a mean error of 4 "7 . From the latter observer, 
23 observations 1836 January 14 to May 17 are available, but for the 
largest part the comparison stars are not contained in modern star 
catalogues. However, after a preliminary re-reduction using star 
positions partly from earlier sources, the mean error of about the 
half of the observations was 5'.'4 so that principially this observational 
serie can be used. For the present investigations, however, it 
appeared senseful to use exclusively the rather precise observations by 
f. w.Bessel (1844) at Kb'nigsberg and F.G.w.struve (1839) at Dorpat 
(two very experienced observers) instead of a lot of much more uncertain 
or only limited re-reduceable observations. An attempt was made to 
re-reduce their observations by means of modern comparison star positions. 
But first, only the half of the used comparison stars are contained in 
modern star catalogues, and secondly, after the reduction and a subsequent 
orbit fit the residuals of these observations have been increased, so 
that obviously the star positions computed back over one century are worse 
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than the positions obtained with special effort in that time by the 
two observers. Thus, finally 27 normal places performed by 
H.westphaien (1847) were used and reduced from the system of the 
36 clock stars to FK^. They received weights of 0.6 each. 

From 1759, 150 observations are published which scatters by about 
0!6 . With completely sufficient accuracy, differences to an ephemeris 
obtained during an earlier investigation (o.A.Rosenberger, 1830) have 
been plotted and the ephemeris corrections for three normals taken from 
this figure. The calculated mean errors are 4" to 8" and thus weights 
of 0.1 in mean are adjoined. 

For the appearence 1682 only measurements of distances of the comet 
to surrounding stars are published instead of equatoreal position values, 
so that here three normals were generated by computing them from orbital 
elements which have been derivated earlier from these measurements by 
o.A.Rosenberger (1830). Because the accuracy of the measured distances 
is of order 1', this orbit and the normals represent them without loss 
of information. The weights have been choosen to be 0.1 . Similar 
like to 1759, the times of the three normal places were choosen under 
consideration of the distribution of the observations. 

For 1607, one normal place was generated by means of ten observations 
reduced by F.t/.Bessei (1804). Two observations each in a and <S and the 
last observation completely have been ignored, and the time of the first 
observation was corrected (T.Kiang, 1972). The corresponding weight 
was calculated to 0,0004, but the perihelion time is determined to 
1607 October 27.5196 +070043 TDB which strongly suggests taking this 
appearence into account in the computations, too. 

For the present apparition, all observations published at the time of 
the individual computations were processed, up to those given on 



IAU-Circular No. 3914 used for the most recent results. Because of 
their uncertainty, the observations of the recovery night made at 
Palomar Mountain, and on 1983 December 31 at Hawaii were not used. 
Although they have lower residual noise, all these observations at 
larger heliocentric distances received one weight unit only, so as 
to avoid forging too many the results due to the errors of the used 
nongravitational force models (cf. chapter 2 and 3 below). 

Besides these standard data, there exist a collection of 663 
individual observations from 1759 on (cf. table 10), which is in use by the 
'International Halley Watch' (IHW) as base for its computations. 
These data correspond, with a few exceptions, to the values originally 
published by the observers, and have been reduced to the equinox 1950 
and into astrometric positions, but without any correction to a common 
reference system. Moreover, accuracy and weight of the different 
observational series are not valued sufficiently. Because of these and 
some other considerations, these data were not prefered about the 
normal places explained above, and have been used only for comparisonal 
purposes. 

2. Some aspects refered to orbit determinations of comets 

With this chapter it is intend to examine some problems in connection 
to orbit computations of comets, with especial regard to the present 
work. 

2. 1 . Nongravitational forces 

Besides the attraction by sun and the planets, the motion of the 
comets is also influenced essentially by the repulsion forces due to 
gas sublimation on the surface of the cometary nucleus. Quantitative 
theories for the forces are given by F.w.Bessel (1836), 
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B.G.Marsden and Z.Sekanina (1968-1972) and especially for COmet 
Halley by H.Rickman and c.Froeschle (1982). The components of the 
nongravitational acceleration can be written in the form 



Ik = g.(r)-A i (t,r) ;; 



(t =1,2,3 for the radial, toroidal and normal direction to the orbital 
motion). The g-(r) should be choice so that they contain the 
dependence of the acceleration components on the heliocentric distance 
completely, and should be normalized so that g.(r = 1 A.U.) = 1. 
Then the A. correspond to the nongravitational acceleration components 
at r = 1 A.U. and only include an explicite (secular) time dependence 
of the nongravitational forces. 

The accurate shape of the g,(r) depends on the direction of the 
sublimation center at the nucleus and thus on the orientation of the 
rotational axis, on the rotational period, and on the thermal inertia 
of the cometary nucleus, as the sublimation center is shifted from the 
subsolar point by an lag angle a(r) in direction of the rotation. This, 
in general, causes an asymmetrical lapse of the g.(r) before and after 
the perihelion transit, if not then the orbital normal vector and the 
rotational axis lie in one plane at perihelion transit. For the 
individual comets, these circumstances in particular are widely unknown, 
and thus only a very rough overall treatment on computation of the 
nongravitational forces is presently possible. 

Bessel assumed constant values of the g, and A i during short time spans. 
Marsden and Sekanina equalized g-(r) with the sublimation rate of water 
ice according to A.H.Delsemme (1982). This implies a rotation axis 
perpendicular to the orbital plane and thus a constant lag angle 
a = arctan (A 2 /^\) and A 3 = 0. However, in general, A 3 4 and a 
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secular variation A. * e~ i£: was considered. Rickman and Froeschle 
i 

carried out model computations for several assumptions about rotational 
period and axis orientation, chemical and physical composition and 
other parameters. From the local surface temperature and the resulting 
sublimation rate and velocity, the magnitude and direction of the 
nongravitational accelerations b-(r) is obtained by integration over 
the whole surface. For seperation of the A., the g.(r) ad hoc 
were equilized to Deisemme's formula, too, so that here the A. depend 
considerably on r. Except for uncertainties in the other assumptions, 
the results depend mainly on the thermal inertia I., of the cometary 
nucleus. This especially is the case for the ratio A 2 (r)/Ai(r) 
whose average value can be determined from positional observations. 
Furthermore, these model computations gave negligible effects into the 
the orbital motion due to A 3 . 

The nongravitational forces produce difficulties on the orbit 
computations due to the following reasons, 

I. The models quoted above satisfy only marginally the entangling 
circumstances. On the theory by Marsden and sekanina, the usually 
more or less unknown rotational parameters remains a priori out of 
account , and a point of large uncertainty on the more explicite models. 
The chemical composition and the accurate surface temperature distribution, 
and correspondingly magnitude and direction of the sublimation are known 
only very approximately. On the models by Rickman and Froeschle, in the 
present state also further essential effects have not been taken into 
account, e.g. a dust layer on the surface, multiple scattering within the 
coma, and inertia! hot sources in the nucleus. Observations from 1984, 
1910 and 1835, especially the careful observations by e.G.w.stmve (1S39), 
show evidence that the widely accepted rotation period of 10 , which has 



been adopted for a part of the model computations, probably is much 
too small. One of the above authors (H.Rickman) kindly communicated 
to the author the results of model computations using 50 rotation 
period, and remarked that further model computations under consideration 
of the explained and further new aspects are in work. 

Besides the global dependence of the nongravitational forces on the 
heliocentric distance, also the essence and magnitude of short- and 
long-term fluctuations of the forces are unknown, but principially we 
have to expect such in connection with the observed optical activities 
like magnitude bursts, jets etc. Considerable activities on comet 
Halley were observed in 1910 and 1835, and these, together with reports 
from earlier apparitions about a tail division, suggest considerable 
activity and essential changes on this comet. On the other hand, the 
nongravitational parameters are comparatively small - the mean acceleration 
is of order 100 m/d 2 - so that fluctuations even of five times the 
averaged forces would need about one month to produce an observable 
position shift of order 1", and this would, for the most part, be canceled 
by continued observations and orbital fits. Because of such fluctuations, 
if they do not become a part of the global models of the forces, we have 
to consider much shorter durations as they produce no positional errors 
of importance (cf. chapter 5.2.1). 

II. Even if the b,(r,t) would be accurately known, difficulties would 
appear in the estimation of the parameters we have to compute, e.g. Ai 
and A 2 according to Marsden and sekanina., or in the necessary correction 
coefficients of the A^r), A 2 (r) by Rickman and Froeschle or similar 
models, which hereafter also are designated with Aj and A2. 

To better explain these circumstances, one might reflect upon table 1. 
It refers particulary to nongravitational forces according to DeJsemme ' s 
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formula, however, the essential conclusions are qualitatively valid 
for any symmetric force model. The nongravitational perturbations 
of the single orbital elements during one revolution can ensue before 
and after the perihelion in the same or in opposite sense. In the latter 
case, the difference in the value on time of perihelion transit at one 
and the same heliocentric distance is the same before and after it, 
so that the nongravitational forces only produce an inequal motion during 
one revolution. In the first case, however, contrarily is produced a 
remaining secular change at each revolution. 

Table 2 gives the corresponding values for comet Halley. As we have 
to expect for any central force b = r with n ^ 1,-2 , the only secular 
perturbation by A x is a perihelion motion. The small change of the 
perihelion time we have to interprete as the duration which the comet 
need to pass this 1'.'4 perihelion shift. A 2 , however, produces an increase 
of q and e on each revolution, corresponding to a delay of 4.2 days each 
revolution. Temporary changes during one revolution caused by both Aj 
and A 2 arise in all elements. The changes due to A x in q and e from perihel to 
aphel, which disappear until the next perihel again, amounts approximately 
to three times of the corresponding changes due to A 2 , which again 
corresponds to the half of the delay of four days until the next perihel. 
Thus, Ai produces a considerable deviation of a few days in the motion of 
the comet near its aphel compared with the unperturbed motion; in case of 
a negative value of A! the comet is too late. The intristic cause is, that 
in large distances the radial nongravitational force component acts nearly 
parallel to the motion of the comet, and thus any positive A x accelerates 
the comet towards its aphel until that is reaches, but after this it 
counteracts its free fall to sun. In smaller heliocentric distances 
(r < 3 A.U.), the temporarly changes in the motion by both Aj and A 2 
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are only poor perceptible and vastly representable by a slightly changed 
value of the excentricity. In practice, this causes a strong correlation 
of values determined for e, A x and A 2 if we have only observations in 
close heliocentric distances. 

During these considerations it was assumed, that the nongravitational 
forces lapse symmetricaly in the ascending and descending part of the 
orbit. If this is not the case, also the perturbations in the elements 
progress asymmetrically, and, especially, a secular change in the 
revolution time produced by A : must be expected. On comet Halley, 
if Ai before the perihel is by 0.01 larger than after it, this effect 
already would amount to +0.58 days. 

Moreover, we might shortly consider the perpendicular force parameter 
A 3 . If the excentricity is not small, the perturbations mainly happens 
close the perihelion and on each revolution in the same direction, so 
that we have to expect secular changes of the orbital plane orientation. 
The temporarly changes due to A 3 are much smaller than these by Ax and A 2 , 
the difference of q and e on the aphel compared with the values on perihel 
only amounts to +1.9-10" 9 and +6.4-10" 9 per A 3 /0.10 , respectively. 
The perturbations in a and i before and after the perihel lapse differently, 
because the nodes are not located symmetrically to the apside line (w=112°). 
The secular perturbations per revolution amounts to Aw = +1?97, Afi = +2!'07 
and Ai = -0V25 per A 3 /0.10 . The change in the arguement of perihel 
is essentially changed by that of the node; in a resting reference frame 
the perihelion moves essentially perpendicularly to the orbital plane, 
corresponding to Ai. 

For the determination of the nongravitational parameters in case of 
certainly advanced g^r), the above considerations allow the following 
conclusions to be drawn. A 2 is well determined with a high degree of 
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accuracy by three or more observed appearences of the comet and the 
perihelion times implied thereby, because a secular change of the 
revolution period is explainable neither by the classical orbital elements 
nor by Ax (assuming symmetric g.), and thus no correlation of A 2 with 
other unknowns occurs. 

For comet Hal ley, by the perihelion times in 1759, 1835 and 1910 
the increase of the revolution period of about four days is determined 
accurately to a few minutes, which corresponds to a relative accuracy 
of 0.1% in Ai. This effect is the most essential nongravitational 
effect in the motion of comet Hal ley (P.H.cowell and A.c.D.Crommelin 

I910, T.Kiang 1972) . 

Principially, Aj can be computed most accurate from the perihelion 
shift between at least two apparitions. On comet Hal ley, however, 
this is not possible with sufficient accuracy before the perihelion 
transit in 1985, because the longitude of the perihel in 1910 is 
determined with a mean error of iO'.'S, that in 1835 by ±1"4, and thus 
presently the uncertainty still is of order of the perihelion shift we 
have to expect for an amount of Ax = 0.10 . Although in practice, 
in case of accurate calculation of the equations of condition this 
effect certainly is taken into account, too, at the present state the 
determination and the results of A x mainly depend on the few recent 
observations in large heliocentric distances. Whilst a seperation 
of Ai and e is always uncertain if observations are available only from 
low heliocentric distances, thanks to these far distanted observations 
this becomes possible because the considerable deviation from mean motion 
in larger distances, as explained above. In case of symmetric g,, the 
analogous effect caused by A 2 do not make trouble because its certain 
knowledge. 
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Because a motion of the orbital plane is not obtainable by another one 
of the unknowns, A 3 is determinable without principial difficulties 
if it is not very small and produces only unobservable effects. However, 
until there remain small residuals due to our insufficient force models, 
especially at large distances and not lapsing along the line of variation 
precisely, on orbit determinations it can easily happen that these 
will be partly compensated by a small change of the orbital orientation 
for the different appearences and thus by a falsified value for A 3 . 
Apart from this, in case of comet Halley the limits of accuracy for A 3 
mainly are set by that of the available observations. 

III. The present situation in practice is, however, that neither 
the true lapse of the g-(r) is known, nor have we enough observational 
information to seperate well all unknowns even in case of an ascertained 
lapse. There are many examples for periodic comets in which, after 
orbit improvements, systematic residuals remain. This strongly permit 
doubts on the adopted force law, even if these observations provide only 
marginal information for improvements only. 

The fact that, besides the rather poor determination of A x especially 
on using observations from 1835 to 1984 only, we have to consider also 
errors of the adopted model for the b.(r) and its assumed symmetry, 
changes the aspects of parameter estimation discussed above essentially. 
First, the observed increase of 4.2 days in the revolution period must 
not necessarily be caused by A 2 , but can also originate partly by Aj . 
Thus, A 2 is no longer very precisely determined and, similary, its 
influence on the observations at large distances what we need for the 
separation of A : . Whilst in case of knowledge of the true force lapse, 
the far observations 1982 - 4 would be very important for seperate Aj and e, 
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in case of considerable erroneously models they become vast valueless, 
because Aj or its values acting at, and also these computed from, 
large and low distances have nothing to do with each other. Because 
our lack of ability to recognize the accurate development of the forces 
until the perihelion transit in 1986, these far observations are presently 
inapplicable for accurate predictions. This is the result of fitting 
a wrong model to these observations, as we then would have to expect a 
corresponding error in the result. 

These essential differences between the previously assumed ideal but 
not present conditions must be admitted. On comet Hal ley we presently 
have the situation that sometimes from observations 1835, 1909 - 11 and 
1982 - 4 , elements inclusive the perihelion time for 1986 have been 
computed and the results published. These observations are quite 
sufficient to seperate the unknowns within a certain force model, so 
that it is not surprising that there remains no systematical residuals. 
However, first this does not permit the conclusion that thus the model used 
ad hoc (Deisenme' s formula) certainly reflects the true force lapse well. 
Secondly, the obtained low mean errors calculated for the unknowns are 
only correct within a certain model and yield of a too high accuracy of 
the results. The formal mean error of the perihelion time in 1986 
appears to ±0.008 days, but the full error of the model enters in the 
calculated perihelion time as explained above, so that its true error can 
be ten times larger. If we use the observations of three or more 
previous apparitions, the perihelion time follows much more precisely 
and rather independently on the adopted force laws and their errors, and 
differences between the far observations permit conclusions about the 
favourable ones of the models instead about the perihelion time in 1986. 
In fact, presently the latter is not better known than before the recovery 
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of the comet in 1982. 

Any quasistatic force model, which assumes the same force lapse on 
each revolution without secular or essential short-term changes, 
produces independently on its explicite form on every revolution the 
same delay in the subsequent perihelion transit time. On comet Halley, 
however, the amount of 4.2 days does not remain constant, but increased 
by about one hour each revolution, as followed significantly by the 
observations back to at least 1531. This cannot be explained due to our 
inability to recognize the force lapse, but suggests a secular increase 
of the nongravitational forces by 1.1% per revolution. Thus, a further 
parameter B has been introduced, in order to describe the time dependence 
by 



A.(t) •= A. (epoch)- (1-B-t) 2) 



for which sign and time unit (10000 days from the epoch) have been 
choosen in accordance to a similar parameter used by Maxsden and sekanina. 
In order to avoid a secular change in the lag angle, as well as in the 
delay of one hour which on the asymmetrical force models can be caused 
by both hi and A 2 , B has been refered to both parameters commonly. 
During application of different models to the observations, the 
results for B changed only minutely and were always determined well. 
There are some physical reasons which let us expect such a change of the 
nongravitational parameters. Due to sublimation, on each appearence the 
cometary nucleus decreases in size. If its constitution, and thus the 
sublimation rate per surface unit element, is assumed to being constant, 
the nongravitational acceleration increases with the ratio of surface 
to mass or volumina, or indirectly proportional to the radius of the nucleus. 
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In the case of a radius of 3 km, the sublimation of a 30 m thick 
layer per revolution, which is necessary to explain the force increasement 
of 1%, appears somewhat too large. However, at least 2Q% or so of the 
observed effect, or about ten minutes in the perihelion time, must be 
expected. Another effect is, that with progressing decay of the comet, 
the lag angle increases. If we assume - for a very rough valuation - 
that A 2 « hi and that the increase of the ratio A 2 /Ai by a factor ten, 
as observed on many old comets, happens within thousand revolutions, there 
would appear an increase of A 2 and of the corresponding revolution time 
increment of IX per revolution. Certainly, under bad circumstances 
this conclusion is invalid (e.g., if the force law is very assymetrical 
and thus the revolution time delay mainly is caused by Aj instead of A 2 , 
we would have to expect an decrease instead), but in any case it is 
obvious that, in general, secular changes of the forces of this amount 
are possible and have to be considered. A further cause can be a 
change of the inclination between the orbital plane and the cometary 
equator, caused either by planetary perturbations of the orbit or by 
precession of the nucleus. The different reasons which allow us to 
expect changes of one hour in the nongravitational revolution time delay 
are a further strong approach to prefer solutions from observations 1607 
to 1984 before such from 1835 to 1984, because rather independent on 
the particularly physical reason, the half of this effect (corresponding 
to a half hour in the perihelion time) still will be originated until 
February 1986 and thus is not contained already in the 1982 - 4 observations. 
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2.2. Displacement of the light center on observations to the 

nucleus of the comet 

If we consider the observational residuals after orbit improvements 
of comets, systematic values during some time intervals become visible 
in some cases and are not explainable in the usual way, e.g. by errors 
in the comparison star positions. Usually this is interpreted as the 
presence of an offset between either the image center or its brightest 
point measured, to the cometary nucleus, which we have to expect because 
the gas emission and points preferably in direction towards to sun. 
However, neither quantitative investigations nor a theory of this effect 
has been presented. The residuals, however, permit the following 
conclusions: a) The light shift approximately points along the radius 
vector to sun. b) Its magnitude is not in simple relationship to the 
heliocentric distance, but obviously depends strongly on the observers 
and on observational circumstances. On very short exposed plates, 
the effect in whole is smaller than on longer exposed ones, and in many 
cases none systematical residuals appear in observation series with very 
short exposure times and stellar images of the comet. But also it 
might be important, if the position was obtained by a densitiometric 
measurement which gives always the brightest point of the coma nearly 
independent on the exposure time, or by visual measurement of the geometric 
center of the image. In the latter case, in addition to the dependence 
on the exposure time, physiological errors along symmetry direction of 
the comet (or correspondingly approximately along the radius vector to sun) 
might possibly occur. Probably there is also a dependence on the spectral 
range of the exposure because the different positions of brightness center. 
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Although in medium-term sight we do not expect that these 
circumstances will be cleared in detail, the question arises 
if this effect cannot be taken into account anyhow during orbit 
computations. In the literature, an offset along the heliocentric 
radius vector is usually assumed, and for its magnitude one adopts 
a radius dependence S(r) = S -s(r), where S is an unknown parameter 
which we have to determine together with the orbital elements using 
different assumptions for s(r). Because of the above considerations, 
however, such an onset probably is rather valueless as the dependence 
on r is obviously much less than that on other facts. 

A possibility to eliminate the effect rather independently in 
assumptions on its magnitude, is a transformation of the equations of 
condition for the orbit improvement from right ascension and declination 
to position angle and appearent angular distance to sun, and the 
exclusive use of the equations for the position angles only. The latter 
ones are not influenced by any radial offsets, whilst the final residuals 
in the angular distances to sun reflect the magnitude of the offset in 
the single observations and perhaps permit conclusions on its dependence 
from observational circumstances. The transformation of the equations 
of conditions to position angle 8 and anqular distance if» can most 
simply be performed by 

(do/ A l U -b \ /da cos s\ 

\ d */ -?p7b?\ b *}{ d& I 3) 



with 



, _ 3a cos 6 _ y cos a - x sin a ., 

3 33 Z? 4 ) 

h _ 36. _ (-x cos a + y sin a) sin 6 + z cos 6 ,-\ 

3S " Ar *> 



(x,y,z,r heliocentric coordinates and distance of the comet), 
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or immediately by numerical elimination of S in the both equations 
for Aa cos 5 and A6 of each observation considered, but in 
regard to the normalisation of the weights. 

By application of this method, the mean errors of the elements 
increased a few times. This is rather unchanged even if both 
positional values are used for some observations in larger 
heliocentric distances. Using observations from two or more 
appearences, the mean errors increased only minutely, most that of 
the perihelion time, and only in case of more than three apparitions 
independently if nongravitational parameters are estimated or not. 
The cause is, that by more than one apparition, the semimajor axis, 
and by the position angles during the single appearences the angular 
motion in the orbital plane and thus the perihelion distance and 
excentricity are determined well, better determinations than these 
being possible only by observations of the true or appearent distance 
to sun. Hence, the elimination of the offset is possible, as the 
concerned distances enter only minutely into the computations in case 
of more than one apparition. However, this is only correct for 
sufficient orbital inclinations, because on the limit sin i ■+ 
the position angles to sun become meaningless and thus on low 
inclinations the uncertainty in e,q and especially T would be 
increased considerably. Because the uncertainty of predicted 
ephemeris places mainly depends on that of the perihalion time, then 
it also can increase several times. In practice, we have to check 
individually if, in the considered case, the increase in position 
uncertainty on ignoring the elongations becomes larger than the 
light offset we have to expect. 
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A remaining question is, how accurate the assumption of a light 
offset along the radial direction is at all. Presently, at least 
based on the available astrometric data of comets, this cannot be 
ascertained better than over the indication quoted above, and thus 
this assumption must be used in order to having a working hypothesis. 
Perhaps in the future the application of different assumptions to 
observations of comet Halley on smaller heliocentric distances will 
give improved verdicts. 

3. Further details on the performance of the computations 

The used observations had been discussed in section 1, and 
here it is intended to explain some further bases and details of 
the computations. 

The computations were carried out by means of a computer programme 
which the author has performed six years ago and which has established 
itself well since without essential modifications. First, initial 
values of the planets are read in and are integrated to the osculation 
epoch choosen for the comet. To integrate the planets, too, has of 
course the disadvantage of much increased computation time, but on 
the other hand it is not necessary to use tapes containing planetary 
coordinates in case of long-term integrations and, moreover, the 
procedure remains flexible in regard to the integration step width 
and to the choosen planetary theory. Afterwards, preliminary values 
for orbital elements of the comet and further unknown parameters, e.g. 
the nongravitational parameters, are read in, and also the observations. 
These are reduced to the FK 5 system so far as presently possible by 
correction of the main term of elliptic aberration, of the equinox and 
precession, and for the proper local errors of the FK 4 system according 

23 



to modern meridian observations of stars. For the nongravitational 
force lapse in the programme different theories are incooperated 
which are coiceable by an input digit (presently those according 

to Marsden and Sekanina / Delsemme and Rickman and Froeschle) and 

which can easily be exchanged for other ones on demand. Values 
of the parameters A : ,6^2.82^3, B 3 , the orbital elements, and 
further ones, or of any linear combination of these, can be either 
considered or estimated. The computation of the coefficients of 
the equations of condition have to be performed rather precisely 
because of the strong correlation between some of the unknowns, as 
explained in chapter 2. A method of great flexibility and accuracy, 
well established also in difficult cases and on a large number of 
unknowns, which calculates the partials by integration simultaneously 
together with that of the celestial bodies, was used for this 
(Landgraf, 1983b). In order to avoid large residuals and to support 
the convergence during the searching of sufficient initial 
values of the unknowns for each force model, in most cases 
the osculation epoch has been placed on 1910. First was integrated 
forewards, then backwards, to compute both the cometary position 
and the equations of condition for all observation times decreased 
for the light time. An integration step size of 0.625 days was found 
to being most favourable for integrations back to 837, and used 
throughout. For solving the conditional equations, the Method of 
Least Squares was adopted. The complete procedure was repeated 
until a certain accuracy or number of iterations was attained. 
The used criterion was the comparison of the remaining error square 
sum guessed during the solution of the normal equations with the 
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value later obtained by the new residuals of the observations. 
Because of the accurate computation of the conditional equations, 
the convergence was \/ery good. Also in cases of rather bad 
starting values, only rarely more than two iterations and one 
computation of the equations was necessary. 

For the planets Mercury to Neptune, initial values obtained 
by the Institute for Theoretical Astronomy at Leningrad (ITA) 
have been Used (G.A.Krasinsky, E.V.Pitjeva, M.L.Sveshnikov and 
E.s.sveshnikova 1982) and refered to the FK 5 system and reciproce 
mass values for Jupiter, Saturn and Uranus of 1047.348, 3498.0 and 
23030, respectively. For comparison purposes, initial values of 
the theory DEI 19 by the Jet Propulsion Laboratory at Pasadena (JPL) 
have been applied. In the first case, standard coordinates of the 
Schwarzschild metric and in the second case isotropic coordinates 
are used, so that the motion of the comet was integrated in, and 
the presented results refered to, the corresponding ones. The 
position of the moon have been computed geocentrically, and the 
perturbations by earth and moon seperately taken into account to all 
other bodies and the earth-moon barycenter. For the transformation 
of the observation times from UT to TDB, values obtained recently by 
comparison of the lunar theory LE200 with observations have been 
kindly communicated by p.K.Seideimann, us Naval observatory, Washington. 

4. Results of the computations 

Since the commencement of the work on comet Hal ley, a large number 
of trials and computations were carried out during which the influence 
of the unknowns on each other and these on the bases have been lighted 
out, and several perceptions for further work have been collected. 
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However, because of the expense, it neither was possible nor 
necessary to always repeat previous computations after some 
additional observations were published. Thus, a few of the 
results represented below are refered not to the most recent 
state of the available observations, but this does not essentially 
defeat the mainly conclusions of these computations. 

4.1, Elements refered to nongravitational forces according to 
the theory by Marsden and Sekanina (Formula by Delsemme) 

The last results obtained during application of the formula 
by Delsemme for the force lapse are given in table 3. The 
elements no. 1 and 3 have already been published in mpc 8665. 

The first two orbits have been computed using observations 
from 1607 to 1984. In contrary to the solutions based on 
1835 - 1984, the observations in 1984 exhibit a systematic residual 
of -l'.'S in right ascension and +0'.'7 in declination in mean. 
However, under consideration of the connections explained in chapter 
2, we must assume that the perihelion time in 1986 is determined 
better by the previous perihelion transits than by far distanced 
observations connected with an ad hoc assumption about the force 
lapse until the perihel, and that any residuals have to traced 
back to insufficiency of the latter. This also is confirmed by 
the fact, that, on adopting other models (cf. chapter 4.2) , 
these residuals decrease by +0'.'4 and -0'.'3 respectively in both 
coordinates, while the corresponding perihelion transit predictions 
based on the same observations differs only within a few 0.001 
without any obviously dependence on the remaining residuals 
(see table 6). These residuals let us conclude that the nongravitational 
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at larger distances decrease much more heavily than according 
the formula by Deisemme (and also than according to the models 
of Rickman and Froeschie) , in agreement with the observed 
activity of comets at large heliocentric distances, and even 
perihelion distances of about 6 A.U. Because the residuals 
try to decrease the result for the perihelion time, the true 
value might be later than computed, and perhaps falls into the 
interval 1986 February 9.51 to 9.55 TDB. We cannot expect that 
this will be cleared accurately before the comet reachs lower 
distances (r < 3 A.U.). Then, due to the lack of improved 
force models, most precise computations of the observations 
in large distances should be either excluded completely, or, 
the two equations of each far observation should be transformed 
to one equation for the unknowns except of the perihelion time, 
by elimination of the latter. Moreover, solutions from observations 
1835 to 1985 ... may then become preferable because of better 
elimination of systematic errors due to long-term variation of the 
nongravitational parameters. 

On consideration of the previous and following results, the first 
orbit in table 1 must be considered as the most favourable result 
on the present state of the observations. 

The second solution corresponds to the first one, with the 
exception that A 3 have been added as a further parameter. For that, 
however, a solution from observations 1835 to 1984 might give a more 
trustworthly result (orbit no. 6). 

The third orbit is based on observations and normal places 
1835 to 1984 only. In this case, the residuals in 1984 are +0!!2 
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in a and -0V2 in s. Alltogether with the observations 1982-3, 
no systematic residual remains, but, as explained above, this 
is no reason to prefer this solution and the following ones to 
solutions 1607 - 1984. This is confirmed again by the backward 
integration results given in table 5 and their comparison with 
the observations. The differences of -0.5, +0.7 and +8 hours in 
1759, 1582 and 1607, respective, due to an average of 8 by other 
unknowns in solutions 1835 - 1984 does not correspond to the 
values of 1 h/rev. 2 we would expect, but neverthless they clearly 
suggest the need of taking the time dependence of the parameters 
for precise predictions into account. After introduction of B 
as an additional parameter, not only the motion 1607 - 1984 can be 
described without significant differences to the observations, but 
also the well determined perihelion time in 837 {see part 4.4 below 
and table 9) is represented accurate to -0.9 days, so that this 
parameter is vindicated and determined well and in good aggreement 
from several earlier apparitions {cf. table 4). 

The fourth orbit, like all further ones, bases on the ihw data 
for 1835 to 1984, and on adoption of the weights like advanced in 
table 10. On the fifth orbit, the planetary initial values have 
been taken from dei 19. It is remarkable that the improved orbits 
computed in connection to the ITA initial values, now in 1835 
exhibit deviation from the observations by up to Q?5 . Like the 
differences between the backward integration of orbits no. 4 and 5 
according to table 5, this is mainly due to the different mass 
values for Uranus (1:23030 on orbit 4, 1:22960 on orbit 5, cf. also 
part 4.3 below). This has to be taken into account in case of 
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employment of the results presented here in connection with other 
planetary initial values. 

The sixth orbit corresponds to the fourth one but with addition 
of A 3 . The latter is similarly well determined like in orbit 
no. 2, but differs in the result. In cases of essential differences 
between the adopted and the true lapse of nongravitational forces, 
the results for the parameters depend very much on the distribution 
of the used observations and on their residuals. Because in the 
solutions 1607 - 1984 the residuals due to insufficient modellation 
of the forces are partly fitted by changes of A 3 , too, but on the other 
hand only in 1835 and 1910 the orbital plane orientation is ascertained 
well enough to permit a determination of A 3 and thus for this parameter 
(in contrary to the computation of A 2 and B) the earlier apparitions 
do not achieve an improvement, the result for A3 from 1835 - 1984 has 
to prefer to that from solution 2. However, the only thing presently 
ascertained is that the absolute value of A 3 is yery probably less 
than 0.1 , in agreement with the expections considering the low values 
of A! and A 2 . 

The seventh and eight solutions have been carried out under 
consideration of a possible light shift in the observations. In the 
first case, in analogy to the literature, a dependence s(r) only on 
the heliocentric distance of the comet was assumed. For this, the 
formula of Delsemme was choosen again (with s(r = 1 A.U.) = 1). 
The result was S = -326 ±75 km , so that the brightest point of 
the coma is displaced towards to the sun, in agreement with our 
expections. The nongravitational parameters and their mean errors 
changed only minutely, in agreement with the fact that we do not 
expect any strong correlation with S . These results are 
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contradictionarly to those obtained at the European Space operations 

Centre at Darmstadt (esoc). An offset of nearly the same amount, 

but in the opposite direction, and an increase of the mean error 

of Aj by 2.5 times (T.A.mrley, 1984) was computed under somewhat different 

suppositions. In case of orbit 8, at heliocentric distances below 

2.0 A.U. only the position angles of the observed positions to sun 

were used. In agreement with the expections on the use of this 

decreased observation matter, the mean error of an unit weight 

decreased slightly [2%), and that of the different elements increased 

slightly, by up to 1.5 times (for Ai and T) . The differences in the 

residuals of the elongations between this solution and solution no. 3 

amounts to only a few O'.'l . This is much below the strong scattering 

of the residuals of the ihw observations, so that they obviously are 

mainly caused by the different sources of the comparison star 

positions used, or by other errors. Because of this, but also 

because of the fact that single observational circumstances like the 

exposure time are in most cases unknown, it is not possible to draw 

any essential conclusions about the light shift from these data. 

In general, during the progress of the investigations it was 
noted that more and more decreasing values for Ai (down to negative 
values) and for the secular increment of the forces, -B, resulted 
due to an increase of the observations used (either by adding earlier 
apparitions, or additional recent observations). A few examples 
for this are given in table 4. 

The only published orbit which is also based on observations 1835 
to 1984 was computed at esoc (T.A.Morley, 1984). In contrary to our 
solutions no. 4ff . , the four recently measured observations of 1911 
mentioned in chapter 1 have not been used, and the weights differed 
slightly (see table 10). The planetary coordinates have been taken 
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from the theory DEI 18, which corresponds to DEI 19 but is refered 
to the FKt+ system at epoch 1972.5 . These earth coordinates, 
inertially at rest, have been combined with cometary observations 
refered to the rotating FKi, system. Because of the somewhat 
different bases, these results are only very limited comparable 
with ours. The difference of 0.005 in the result for the perihelion 
time to our solutions 4 to 8 corresponds to a position shift of 0V1 
in the 1982 - 4 observations, and thus is within the limits to be 
expected because of the corrections to the FK 5 system on the latter 
solutions. The other elements agree with ours within their noise, and 
also the mean error and the residuals of the observations are in good 
agreement. 

4.2. Elements refered to nongravitational forces according to 
the models by Rickman and Froeschle 

The results obtained on application of the models by Rickman and 

Froeschle (1982 and private communication) are given in table 6; 

A x and A 2 are the necessary correction factors of the model values. 
The results for A 2 correspond to the ratio of the value yielded by 
Deisonme's formula to that according to the applied models at 
r = 0.6 ... 0.8 A.U., so that this interval points out the averaged 
value of the forces. 

All the solutions 1835 - 1984 fit equally well the observations 
at large heliocentric distances. The representation of the earlier 
appearences differs largely. Although the representation is clearly 
better than by the Delsemme formula (no. 3 in table 5), it is not 
possible to prefer some of the models from these results. Considering 
the results for A], and A 2 , the models of higher thermal inertia, which 
also gave the lowest mean residual of the observations, appears to be 
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most favourable. 

By introduction of B and using observations from 1607 to 1984, 
again it was possible to represent all observations well. The good 
agreement of the well determined result for B confirm it's justification 
again. The mean residuals of the observations are only a little better 
than on application of Delsemme's formula, but the representation 
of the far distanced observations is significantly better, especially 
on the models of higher thermal inertia. Because the results for Ai 
have been very low during the application of the rather different models, 
we must conclude that the radial force component is truely negligible. 
Furthermore, on use of these models Aj and B again decrease with 
increasing number of apparitions and observations used. For example, 
using the first one of these models, from observations 1759 - 1983 
resulted A t = +0.43, from 1607 - 1983 A x = +0.12 . 

Altogether, we can conclude that the secular increase of the 
nongravitational forces of comet Halley is well determined, that Aj 
is nearly zero and thus the position errors at far distances are 
mainly caused by an essentially asymmetrical lapse of A 2 instead of 
any Aj , that the nongravitational forces reach to larger distances 
than represented by all available models, and that the mckman - 
Froeschie models are clearly favourable to the Delsemme formula, but 
that also these models still are far away from a representation of 
the true force lapse. For more detailed conclusions, however, further 
observations and models are urgently necessary. 
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4.3. The influence of the masses of Uranus and Neptune 

During the computations it have been noted that, besides other 
influences, the results rather depended on changes of the adopted 
mass values for Uranus and Neptune within the limits of their 
uncertainties. This is caused by the fact that these planets 
have approximately the same and double revolution time like the 
comet, respectively, and that on each or every second one of the 
last stayings of the comet on the far parts of its orbit, these 
planets have been in similar heliocentric direction. 

Table 7 gives the changes in earlier perihelion times corresponding 
to solutions 1759 - 1983, due to variation of the reciproce mass values 
by +50 units and due to variation of B 2 , respectively. At the time 
of performation of these computations the latter have been used, 
and contrarily to the later used B it is only refered to A 2 ; however, 
because of the poor determination of Ai both values coincide. 

Here would be the wrong place to carry through a discussion of 
the most probable mass values for Uranus and Neptune and their accuracy. 
Under consideration of the results of the different determinations 
(a review, for example, is given by L.Ballani, 1981), however, we can 
say that the uncertainty that is to be expected in the mass values for 
Uranus and Neptune, and correspondingly those in the computed perihelion 
time, are approximately one and two times of the range of table 7, 
respectively. A variation of the mass of Uranus within acceptable 
limits can make amends only for a small part of the results of B 2 or B, 
but the vast correlation between both parameters has the practical 
advantage that errors in the adopted mass for Uranus are compensated 
by inclusion of B as a further unknown and use four or more apparitions. 
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From the observations 1607 - 1984 in connection with the perihelion 
time in 837 (see part 4.4), a good separation of all parameters 
considered in table 7 is possible, and the reciproce mass of Uranus 
is presently determined accurately to ±40 units. For reliable results, 
however, it is better to wait until improved force models and observations 
close the present perihelion are available. 

4.4. The perihelion transit in 837 and the long-term motion of 
the comet 

To both check the different force models and investigate the 
long-term motion of the comet, a well determined perihelion time 
of a much earlier apparition would be of very large value. It was 
noted that for this the apparition in 837 can be used. 

The comet passed the earth on 837 April 10.63 TDB at only 0.0325 A.U. 
distance and, because of the differential perturbations, a variation 
of the accurate time of the encounter by only 0.1 day would produce 
differences of some days in the previous perihelion times. In table 8 
the results of backward computations until 141 are compared with the 
observed perihelion times (T.Kiang 1972, i.Hasegawa 1979), starting 
with the assumptions T = 837 February 28.40 and 28.44, respectively. 
All other elements are taken from a backward computation of an earlier 
1607 - 1983 solution (see table 9). Corresponding to the observed 
perihelion times in 607, 530, 374 and 141, the perihelion time in 837 
has been between 837 February 28.43 to 28.48, and the good agreement 
suggests that the influence of possible inequalities in the motion of 
the comet to the result can be only very small. By exclusive use of 
the observed perihelion time in 141, which formaly gives the most accurate 
result and, moreover, must be prefered as base of a continued backward 
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integration, we get 837 February 28.427 TDB, or 28.424 TDB at epoch 
837 March 10, with an accuracy which is very probably better than 
±0.010 days. Because, in particular, not the perihelion time, but 
the circumstances of the close encounter to the earth are determined 
by the previous motion, the above result has to be corrected by 
AT 837 = -0.322- (ir-304?140) , where ir=a-o> refers to the epoch 
837 March 10 and to the FK 5 equinox at B1950 (w.Landgraf, 1983c) . 
This result for the perihelion transit in 837 is slightly later 
than that obtained by the observations which have been made in that 
time (i.Hasegawa, 1979). Backward integrations of several solutions 
1607 to 1983 - 4 gave 837 February 27.1 - 27.6 TDB and thus are 
approximate one day too early. 

Starting from this improved perihelion time for 837, subsequently 
the motion of the comet was computed back to 2300 B.C. The results 
are given in table 9. As is to be expected, the perihelion times back 
to 141 are satisfied completely, and also the well observed appearence 
in 12 B.C. with an difference of four days. The earliest ascertained 
apparition is that of 466 B.C. In Greece, the comet was observed 

th 

in the second year of the 78 olympiad in western direction and 
has been described among others by pline and Aristoteies (s.Lubienietzky, 
1668, A.G.Pingne 1783, a. a. Barret 1978). In China it was observed on the 
second reign year of emporer Ting Wang (p.y.ho, 1962). This corresponds 
to the time spans July 467 to June 466 B.C. and February 467 to January 
466 B.C., respective, so that the comet must have been observed 
between July -466 and January -465. This is in agreement with our 
computed perihelion time, because, according to this, the comet would 
have been observable during the winter in western direction before 
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earning into conjunction with the sun. The comet was also possibly 
observed in 618 B.C., but the corresponding report and its date are 
rather poor (A.cpingre, 1783). For the earlier calculated perihelion 
times, no corresponding reports of an observed comet was found. 
Reliable evidence about the accuracy of the back computations are 
first possible after a repetition using an improved theory for the 
nongravitational forces. Because the comet often very closely 
(to a few 0.01 A.U.) encountered the earth, the back computations 
already in 466 B.C. are possibly so uncertain that a close encounter 
could have occured on this or a previously apparition so that the 
earlier motion happened rather differently than according to 
our integration. 
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Table 1 — Perturbations in the orbital elements by nongravitational 
forces 

The perturbations h(per-orig) of the elements in nearly parabolic 
orbits from great heliocentric distances until the perihel by 
nongravitational parameters A. and A 2 of the model by Delsemme 
are presented. If a (+) or (-} is indicated, the perturbations 
h(fut-per) from the perihelion until the following aphel are of 
the same or different sign, respectively. In the first case, 
a secular effect twice the given amount originates on each perihel, 
in the second case only a temporary perturbation. The proper part of 
the perturbations happens in the part with heliocentric distance 
below 4 AU. 
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Table 2 — Perturbations in the elements of comet Halley by 

nongravitational forces according to Delsemme's formula 

Part a): Secular changes of the elements &(per + -per) between two revolutions. 
Part b): Changes in the elements from perihel until subsequent aphel, 
b(fut-per). The perturbations of the elements are performed mainly below 
3 AU heliocentric distance, whilst the deviation Av of the true anomaly 
from the unperturbed motion, caused by Ae and Aq accumulate until the aphel. 
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Nr. 



T(TDB) 



1936 Feb 

1935 Feb 
1935 Feb 
1935 Feb 

1935 Feb 

1936 Feb 
1985 Feb 
19S6 Feb 



9.50752 
9.50235 
9.44912 
9.45034 



45032 
45103 
44S43 
45394 



0.5S7104S 
0.5371049 
0.5371014 
0.5371033 
0.5371030 
0.5371033 
0.5371041 
0.5371061 



0.96727910 
0.95727841 
0.95727328 
0.96727517 
0.96727592 
0.96727535 
0.96727510 
0.95727545 



111.84718 
111.84768 
111.84651 
111.84590 
111.84703 
111.34654 
111.84635 
111.34720 
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58.14427 
53.14386 
58.14417 
58.14424 
58.14378 
58.14418 
58.14431 



162.23917 
162.23928 
162.23940 
162.23933 
162.23932 
162.23928 
162.23933 
162.23929 



ESOC 1936 Feb 9.44367 0.5371022 0.9672750 111.84690 58.14414 162.23938 



Nr. 



1 


-0.0133 +0.0092 


0.015964 +0.000016 


2 


+0.0192 "0.0100 


0.015935 +0.000016 


3 


+0.1232 +0.0205 


0.015516 +0.C00024 


4 


+0.0763 +0.0203 


0.015519 +0.000017 


5 


+0.0531 +0.0204 


0.015504 +0.000017 


6 


+0.0728 +0.0205 


0.015521 +0.000017 


7 


+0.0763 +0.0204 


0.015512 +0.000017 


8 


+0.0439 +0.0301 


0.015549 +0.000021 


OC 


+0.080 


0.0155 



-0.00495 +0.00018 
-0.00472 +0.00019 



A 3 =+0.0556 +0.0159 



A.=-0.0396 +0.0143 
Sq=-325 +75 km 



Nr. 



arc 



Observations 
no. a 



1607 - 1984 
1607 - 1934 
1335 - 1934 
1335 - 1934 
1335 - 1934 
1335 - 1934 
1335 - 1934 
1835 - 1934 



ESOC 1335 - 1984 



91 

91 

84 

662 

652 

662 

652 

(662) 

653 



1'.'43 
1.43 
0.94 
2.15 
2.15 
2.15 
2.15 
2.11 

2.20 



1'.'25 
1.25 
0.92 
1.02 
1.02 
1.02 
1.02 
1.00 

1.07 



FK5 
FK5 
FK5 
FK5 
FK5 
FK5 
FK5 
FK5 

(FK4) 



ITA 
ITA 
ITA 
ITA 
DE119 
ITA 
ITA 
ITA 

DEI IS 



Epoch 1936 Feb 19.0 TDB, Equinox B1950. 



Table 3 — Orbital elements of comet Bailey with nongravitational forces 
according to Delseame's formula 

K reference coordinates for the elements: N newtonean, S standard- and 

I isotropic coordinates of the Schwarzschild metric 
S reference system 
P adopted initial values for the planets: ITA from ITA Leningrad (Mercury 

to Neptune) ,VlEl I8.DE1 19 from JPL Pasadena (Mercury to Pluto) 
o,U root mean square and unit weight residual of the used observations (cf .table 10) 

See also table 16. 
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Observations 


T(TDB) 


A l 


A 2 


1835 - 1934 


1986 Feb 9.449 


+0.12 


+0.01552 


1607 - 1934 


9.508 


-0.01 


+0.01595 


1607 - 1933 


9.529 


+0.01 


+0.01609 


1759 - 1933 


9.535 


+0.03 


+0.Q1610 


1759 - 1932 


9.474 


+0.13 


+0.01551 


1682 - 1532 


9.549 


+0.10 


+0.01619 





-0.0050 
-0.0059 
(-0.0059 assumed) 


-0.0067 



Epoch 1936 Feb 19.0 TDB 



Table 4 — The dependence of the results for A ,A ,B and T on the 

employed observation matter 
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Epoch (TDB) 



obs 



1986 Feb 19 1986 Feb 9 9.50762 9.44912 9.45034 9.45082 

1910 May 9 1910 Apr 20 20.17871 20.17859 20.17849 20.17852 

1835 Nov 18 1835 Nov 16 16.43953 16.43961 16.43953 16.43956 

1759 Mar 21 1759 Mar 13.0628 +0.0012 13.05932 13,03703 13.04513 13.04255 

1682 Aug 31 1682 Aug 15.2806 +0.0022 15.28158 15.30807 15.32121 15.30713 

1607 Nov 13 1607 Oct 27.5196 +0.0043 27.51776 27.89802 27.91202 27.75982 



Table 5 — Comparison of the perihel times of the single apparitions 

The given no. corresponds to those in table 3. In contradiction to the 
solution 1 front observations 1607 - 1984, the other orbits from observations 
1335 - 1984 exhibit significant differences to the previously perihel times. 
The large differences for 1607 on orbit 5 compared with orbit 3 and 4 are 
mainly caused by the differently adopted mass values for Uranus (cf. table 7). 



Table o — Orbital elements of comet Halley with nongravitational forces 
according to the models by Rickman-Froeschie 

All represented orbits were computed from the same observation matter 

1607 - 1984 and the same further bases like solution 1 in table 3. 

ia and A5 are the difference in the representation of the observations from 

1934 ro that solution. Orbit 1-3 refers to models assuming P =10" 

ar.d I =130,500 and 1000, orbit 4 and 5 to P =50^ and I =130 X and 100O. 
th rot th 
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Table 6 



No. 



a) 1607 - 1984 



T(TDB) 



1 


1986 Feb 9.50517 


0.5871077 


0.96727826 


111*84716 


58114373 


162°. 23921 


2 


1936 Feb 9,50390 


0.5871061 


0.96727842 


111.84716 


58.14370 


162.23920 


3 


1986 Feb 9.50458 


0.5871054 


0.96727858 


111.84716 


58.14367 


162.23919 


4 


1986 Feb 9.50556 


0.5871085 


0.96727780 


111.84716 


58.14374 


162.23922 


5 


1986 Feb 9.50304 


0.5871062 


0.96727839 


111.84716 


58.14371 


162.23921 



No. 



Al 



A 2 



Aa 



AS 



obs . 



1 -0.020 ±0.061 


+3.415 ±0.002 -0.00395 ±0.00019 


+0'.'37 


-0'.'28 


l 1 


126 


91 


2 -0.032 ±0.060 


+1.585 +0.002 -0.00400 ±0.00019 


+0.30 


-0.20 


1 


.25 


91 


3 -0.049 ±0.061 


+1.088 ±0.001 -0.00426 ±0.00018 


+0.19 


-0.12 


1 


.27 


91 


4 +0.036 ±0.063 


+6.445 ±0.006 -0.00390 ±0.00019 


+0.44 


-0.33 


1 


.27 


91 


5 -0.026 ±0.062 


+1.802 ±0.002 -0.00391 ±0.00019 


+0.35 


-0.24 


1 


.24 


91 


Epoch T 


1 2 


3 




4 


5 






observed 


19db 02 19 1986 Feb 


9.50517 9.50i9Q 


9.50458 


9, 


.50556 


9.50304 








1910 05 09 1910 Apr 


20.17869 20.17369 


20.17870 


20. 


.17868 


20.17869 








1835 11 18 1835 Nov 


16.43952 16.43953 


16.43953 


16 


,43952 


16.43953 








1759 03 21 1759 Mar 


13.05989 13.05957 


13.05931 


13 


.06015 


13.05965 




13. 


.0628 ±0.0012 


1682 08 31 1682 Sep 


15.27975 15.28001 


15.27998 


15. 


.28009 


15.28063 




15. 


,2806 ±0.0022 


1607 11 13 1607 Nov 


27.51840 27.52749 


27.51359 


27 


.52247 


27.52035 




27. 


5196 ±0.0043 


b) 1835 - 1984 



















No. 



T(TDB) 



1 


1986 Feb 9.44889 


0.5871029 


0.96727146 


Ul°84632 


58° 14387 


162°23942 


2 


1986 Feb 9.45022 


0.5871021 


0.96727322 


111.84644 


58.14387 


162.23942 


3 


1986 Feb 9.45092 


0.5871017 


0.96727419 


111.84651 


58.14387 


162.23942 


4 


1986 Feb 9.44817 


0.5871034 


0.96727012 


111.84624 


58.14387 


162.23942 


5 


1986 Feb 9.45024 


0.5871021 


0.96727318 


111.84644 


58.14337 


162.23942 



No. 



Al 



A 2 



obs . 



1 


+1.095 ±0.136 


+3.331 ±0.005 


0'.'91 


84 


2 


+0.934 ±0.132 


+1.545 ±0.002 


0.90 


84 


3 


+0.854 ±0.131 


+1.060 ±0.002 


0.89 


84 


4 


+1.251 ±0.140 


+6.300 ±0.009 


0.91 


84 


5 


+0.968 ±0.136 


+1.777 ±0.002 


0.89 


84 



Epoch 

1986 02 19 
1910 05 09 
1835 11 
1759 03 
1682 08 31 
1607 11 13 



13 
21 



1986 Feb 

1910 Apr 

1835 Nov 

1759 Mrz 

1682 Sep 

1607 Nov 



1 

9.44839 
20.17857 
16.43962 
13.02589 
15.27062 
27.71973 



observed 



9.45022 
20.17859 
16.43960 
13.03025 
15.28074 
27.75112 



9.45092 
20.17861 
15.43959 
13.03342 
15.29092 
27.79830 



9.44817 
20.17855 
16.43963 
13.02377 
15.26572 
27.70633 



9.45024 
20.17859 
16.43960 
13.03026 
15.28004 
27.74671 



13.0628 
15.2805 



±0.0012 
±0.0022 



27.5195 ±0.0043 
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Epoch 


-i- +50 


^=+50 


B =+0.0005 


1986 


-0^00715 


-0^00029 


-0^00530 


1682 


+0.00545 


+0.00455 


+0.00588 


1607 


+0.04465 


+0.03566 


+0.04428 


837 


+0.174 


-0.022 


+0.316 



Zaile 7 — Dependence of the perihel times in 1986, 1682, 1607 and 837 
predicted by solutions 1759 - 1983, under variation of the 
mass values for Uranus and Neptune and of B 

Jiven are the changes in perihel time in days after variation of the reciproce 
mass values for Uranus and Neptune by *-50 units, and of B 2 by +0.0005 . 
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I II T , 

obs 

837 Feb 28.40 837 Feb 23.44 

760 May 20.93 760 May 20.47 

584 Oct 2.36 684 Sep 30.98 

607 Mar 14.13 607 Mar 13.29 607 Mar 12.5 +1.5 

530 Sep 25.15 530 Sep 25.90 530 Sep 25.5 +1.5 

451 Jun 26.04 451 Jun 27.88 

374 Feb 14.08 374 Feb 15.93 374 Feb 16 +1.5 

295 Apr 21.98 295 Apr 19.76 

218 May 23.20 218 May 14.16 

141 Mar 31.03 141 Mar 16.27 141 Mar 21 .1 +1.5 



Table 3 — The representation of the apparitions back to 141 on variation 
of the perihelion time in 837 
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a) 



Epoch (TDB) 



T(TDB) 



2284 06 16 
2209 02 10 
2134 03 15 
2061 08 04 
1986 02 19 
1910 05 09 
1335 11 13 
1759 03 21 
1682 08 31 
1607 11 13 
1531 08 14 
1456 06 28 
1378 11 05 
1301 11 09 
1222 10 15 
1145 04 02 
1066 03 08 
989 08 19 
912 07 14 
837 03 10 



2284 06 
2209 02 
2134 03 
2061 07 
1986 02 
1910 04 
1835 
1759 
1682 09 
1607 10 
1531 
1456 
1378 
1301 



11 

03 



1222 09 

1145 04 

1066 03 

989 09 

912 07 

837 02 



06.84231 
05.49280 
28.66037 
28.86064 
09.50762 
20.17871 
16.43953 
13.05932 
15.28158 
27.51776 
26.25S52 
09.49917 
10.62287 
25.18611 
23.55316 
18.11662 
20.06530 
04.09034 
17.00435 
27.53654 



0.59211254 
0.59018100 
0.59322261 
0.59278730 
0.58710485 
0.58720991 
0.58656496 
0.58446927 
0.58262489 
0.58364892 
0.58122493 
0.57974055 
0.57627478 
0.57277443 
0.57429114 
0.57489917 
0.57458503 
0.58201659 
0.53027612 
0.58245958 



0.96681684 
0.96702148 
0.96664479 
0.96657957 
0.96727910 
0.96730558 
0.96739853 
0.96769098 
0.96793295 
0.96750342 
0.96776025 
0.96800529 
0.96838164 
0.96893884 
0.96885062 
0.96378998 
0.96887294 
0.96789357 
0.96807602 
0.96779533 



115.03325 

114.70857 
113.98893 
112.03456 
111.84718 
111.71809 
110.68479 
110.68921 
109.20287 
107.52900 
106.95515 
105.81548 
105.28010 
104.48707 
103.83558 
103.69886 
102.46544 
101.47334 
100.77435 
100.09855 



61^98002 

61.50399 
60.59178 
58.67675 
58.14364 
57.84577 
56.80149 
56.52871 
54.85078 
53.05286 
52.34172 
51.15294 
50.31277 
49.44723 
48.60076 
48.36140 
46.92993 
45.86508 
44.95956 
44.24273 



161=75497 
161.77262 
161.74586 
161.96180 
162.23917 
162.21550 
162.25562 
162.36927 
162.26146 
162.89776 
162.90943 
162.88167 
163.10367 
163.06644 
163.18236 
163.21368 
163.10153 
163.38893 
163.29993 
163.43779 



A x -0.0138 



A 2 +0.015964 B 2 -0.00495 (Epoch 1986 02 19) 
Ecliptic and Equinox B1950. 
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b) 



Epoch (TDB) 


t(t: 


IB) 


q 


e 


(ii 


R 


i 


837 03 


10 


837 02 


28.4241 


0.5824596 


0.9677998 


100910398 


44?24514 


163?43737 


760 02 03 


760 05 


20.615 


0.581957 


0.967864 


100.0022 


44.0035 


163.4330 


684 09 


29 


684 10 01.430 


0.579734 


0.968153 


99.1603 


43.1252 


163.4084 


607 03 


18 


607 03 


13.571 


0.581007 


0.968041 


98.8116 


42.5888 


163.4666 


530 07 


20 


530 09 


25.625 


0.575779 


0.968706 


97.5970 


41.3071 


163.3351 


451 04 


06 


451 06 


27.230 


0.573927 


0.968916 


97.0499 


40.5498 


163.4692 


374 03 


01 


374 02 


15.292 


0.577417 


0.963581 


96.5341 


39.9211 


163.5323 


295 04 


25 


295 04 


20.632 


0.576120 


0.968752 


95.2667 


38.4563 


163.3577 


218 03 


20 


218 05 


17.709 


0.581677 


0.967964 


94.1737 


37.2550 


163.5640 


141 02 


12 


141 03 


21.076 


0.583379 


0.967841 


93.7197 


36.5671 


163.4265 


66 03 


18 


66 01 


21.896 


0.585360 


0.967556 


92.6798 


35.4781 


163.5670 


-11 12 


27 


-11 10 


05.995 


0.587443 


0.967381 


92.5830 


35.2523 


163.5790 


-86 08 


23 


-86 08 


03.536 


0.586875 


0.967670 


90.8193 


33.3840 


163.3324 


-163 12 


25 


-163 10 


30.105 


0.584770 


0.967668 


89.1685 


31.4448 


163.6947 


-239 04 


28 


-239 04 


16.516 


0.585340 


0.967653 


88.0191 


30.0449 


163.4238 


-314 06 


01 


-314 05 


15.216 


0.587653 


0.967311 


86.7582 


28.7325 


163.5810 


-390 03 


12 


-390 04 


28.980 


0.588377 


0.967269 


86.6838 


28.5122 


163.5793 


-455 04 


15 


-465 04 


11.147 


0.588596 


0.967407 


85.0953 


26.8145 


163.2730 


-541 12 


10 


-541 12 


17.112 


0.590201 


0.967168 


83.5615 


25.1016 


163.5488 


-617 09 


20 


-617 09 


19.971 


0.590725 


0.967172 


82.9646 


24.3105 


163.4449 


-692 10 


23 


-691 01 


03.037 


0.590124 


0.967354 


81.5788 


22.8803 


163.3756 


-768 02 


25 


-768 02 


02.102 


0.537300 


0.967652 


81.0201 


22.0952 


163.4349 


-845 06 


29 


-845 05 


20.163 


0.533601 


0.968270 


79.9331 


20.9705 


163.1551 


-923 01 


29 


-923 02 


21.150 


0.584879 


0.968047 


78.7296 


19.4916 


163.3239 


•1001 12 


25 


-1001 10 


12.603 


0.580168 


0.968686 


77.9033 


18.5697 


163.1327 


-1081 10 


26 


-1081 12 


22.048 


0.580099 


0.968692 


77.0967 


17.5479 


163.0869 


•1158 04 


13 


-1158 05 


15.891 


0.586248 


0.967932 


76.6012 


16.9607 


163.1110 


•1236 04 


22 


-1236 04 


08.118 


0.582026 


0.968710 


74.6496 


14.8214 


162.5231 


-1315 06 


16 


-1315 04 


26.872 


0.585688 


0.968162 


72.9303 


12.9420 


162.6629 


•1393 12 


03 


-1393 10 


11.979 


0.585227 


0.968372 


72.5614 


12.4545 


162.5847 


■1472 03 


12 


-1472 03 


03.117 


0.586410 


0.968301 


71.7378 


11.4107 


162.3520 


■1550 08 


29 


-1550 OS 


23.932 


0.592668 


0.967509 


70.8977 


10.4604 


162.3886 


■1627 02 


14 


-1627 01 


07.750 


0.593452 


0.967747 


69.9578 


9.2762 


161.8652 


•1705 08 


03 


-1705 OS 


23.832 


0.598066 


0.957097 


68.8350 


8.0064 


161.9380 


■1782 12 


05 


-1782 11 


03.000 


0.602695 


0.966398 


67.6358 


6.5454 


161.3354 


-1358 09 


15 


-1853 10 


21.969 


0.603160 


0.966699 


66.5011 


5.2945 


161.4205 


-1935 03 


10 


-1935 10 


17.220 


0.604625 


0.966532 


66.3097 


4.8877 


161.3508 


-2009 02 


20 


-2009 01 


06.670 


0.605355 


0.966513 


65.5906 


4.1540 


161.2762 


-2086 01 


15 


-20S6 03 


27.721 


0.605553 


0.966435 


65.2602 


3.5731 


161.2362 


-2162 09 


11 


-2162 10 


31.405 


0.607084 


0.966682 


64.3009 


2.6263 


160.7396 


-2238 01 


13 


-2239 12 


30.106 


0.607594 


0.966544 


62.8200 


0.9533 


160.6292 


-2316 12 


08 


-2316 12 


10.342 


0.609346 


0.966774 


61.8347 


359.7110 


159.9643 






A. +0.055; 


' A 2 +0.01200 B 2 -0.00653 


(Epoch 837 


03 10} 





Ecliptic and Equinox B1950. 



Table 9 — The notion of the comet during 2300 BC to 2300 AD 

The results for the tine span 837 to 2300 AD (part a) base on orbit 1 on 
table 3. The results for 2300 BC until 837 AD bass on a back integration 
of a solution 1607 - 19S3 and correction of the perihel time in 337 (part b) . 
Because the insufficiency of the force model, the uncertainty in the predicted 
perihel tines for 2061 and 2209 are C^l and l d , respective. 
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0.C4 0.2 0.4 1.0 





( 5V5 2VO 


!'.'5 


IVO ) 


a 


u 


<ep> 

r a 


< ep > 5 


1835 


2 173 - 


- 


- 


41'44 


0'.'39 


+1'.*17 


-O'.'Sl 


1910 


1 - 409 


31 


12 


2.20 


1.07 


+0.29 


+0.17 


1986 


6 - 


- 


37 


0.57 


0,57 


-0.09 


-0.09 



1835 - 1984 9 173 409 31 49 2.15 1.02 +0.23 +0.07 



Table 10 — Shore error analysis of the IHSf data for 1835 to 1984 

For the apparitions 1835,1910,1936 and all observations together, is 
given the distribution over different weight classes (in parentheses 
enclosed the mean error a according to IHH and used by SSOC on the elements 
given in table 3), the root mean square and a posteriori unit weight 
mean error, c^-Ee^/Ip and u 2= EE Z /(n-u) , and the average residuals 
in a and 5 . 
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5. The ephemeris uncertainty for comet Halley 
in March 1986 

With regard to the space missions to comet Halley in 
March 1986, it is not without interest to know the position 
of the comet and its uncertainty. Below the uncertainty 
we have to expect in the position determinations for 
13 March 1986 is considered because at this time the 
Giotto spacecraft of esa will be targetted with high 
precision to flyby the nucleus of comet Halley at a 
distance of 500 km. 

5.1. Mean error (variance) of the estimated position 
values 

The mean errors of unknown parameters (e.g., T,q,e,u,fi,i , 
Ai.Aj) to be determined from observations (e.n. a, <5 ) , 
as well as the mean errors of functions of the unknown 
parameters {e.g. the ephemeris place at the flyby or 
the miss vector), depend on the accuracy of the observations 
and the functional dependence of the parameters from the 
observations. In particular, the mean errors depend 
on the assumed model, the equations of conditions in it, 
and on the number, distribution and accuracy of the 
observations, but not on their accurate values or residuals. 
Therefore, the uncertainty in the position predictions in 
March 1986 can be predetermined, but the results depend 
essentially on the assumptions about the observations until 
March 1986. 
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In table 11 are given some different assumptions 
for the observations until 1986, which are the bases 
for the subsequent error estimations. Six different 
cases assuming different sums of weights are considered, 
from a large number of observations to very bad expections. 
An unit weight corresponds to one observation with a mean 
error of tl'.'O, and given is the assumed sum of weights 
collected for time spans up to one month and placed to 
the date of greatest elongation between the comet and the 
moon . 

Case no. 1 represents most favourable expections. It was 
assumed that the lapse of the nongravi tational forces is 
known with sufficient accuracy, so that the systematic error 
in representation of the observations in large heliocentric 
distances is much below the mean error of its whole 
(approximately 0'.'2), and thus these observations can be used 
for the orbit determinations. Furthermore it was assumed, 
that no light shift exists between the nucleus of the comet 
and the neasured positions, so that also the observations in 
small heliocentric distances can be used. The total weight 
was choosen to approximately 1200 weight units, based on the 
number of observations obtained during the last apparition of 
the comet. 

In case 2, the same assumptions on the nongravi tational 
modelling and the light offset were made, but fewer observations 
were assumed. In addition to the observations already existing 
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until February 1984, a weight sum of 260 till end 1985 
and of additional 6 for the beginning of March 1986 
was assumed. 

In case 3 it was assumed that the modelling of the forces 
is too uncertain for the observations at large distances can 
be used. The limit was set to r = 2.8 A.U. Similar to 
case 2, a low number of observations was assumed, with the 
exception of about twenty weight units for the beginning 
of March 1936. 

The three following cases should investigate the situation 
if one have to consider a radial light offset at low helio- 
centric distances. Similar to a stiff assumption of a certain 
lapse for the nongravitational forces, by assuming a certain 
dependence on heliocentric distance, e.g. Ar(r) = S„'s(r) 
and solving for S„ as an additional parameter, the resulting 
increase in position uncertainty is reflected only insufficiently, 
The results presented subsequently refer to the assumption, that 
nothing is known about the magnitude of the light shift, and 
thus only position angles of the corresponding observations 
will be used (cf. chapter 2.2). The position error is probably 
even more underestimated also by this method, because of the 
possible additional presence of a tangential light shift. 
However, nothing better is presently obtainable, because, 
depending on whether the error of this assumption is erratic 
or systematic, either it must be taken into account by 
increasing the mean observational error or it cannot be 
considered within the postulates of error computations. 
Subsequently, for r < 2.0 A.U. only the weights for the 
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position angles to the sun have been used. Case 4 corresponds 
to the same other assumptions as in case 2. Case 5 
corresponds to case 3 but with the exception of 30 weight 
units in total for March 1986. Case 6 finally assumes only 
6 weight units in March 1986, and shall represent the worst case. 

In addition to these assumed observations in the present 
apparition, the observations of the previous ones 1835 - 1910 
or 1607 - 1910 were added, respectively. On the solutions 
1835 - 1986, the unknowns which had to be estimated are the six 
orbital elements as well as Ai and A 2 ; on solutions 1607 - 1986 
also B (for details, see chapter 3 and 4). In order not to 
repeat in each case the computation of their equations of 
condition, all the observations of the earlier apparitions in 
all cases were used (also on r > 2.8 A.U., and for r < 2.0 A.U. 
all elongations from sun). The thereby slightly decreased 
mean error in the earlier perihelion times and other elements 
have only minute influence on that of the present apparition. 

The results of the error estimations are collected in table 2. 
Under very good conditions, one can expect mean errors of below 
50 km in the impact plane of Giotto and of 120 km along its 
flight direction. In general one has to expect greater 
uncertainties, but it was noted that any omitting of observations 
in r > 2.8 A.U. can be easily compensated by increasinn 
observational effort in March 1986. Also, in the most 
favourable cases, the calculated uncertainty for right ascensions 
and declinations of about O'.'l is above the correlating errors 
one might expect in observations of shorter time spans (assuming 
the use of the new Halley comparison star positions), so that 
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the assumed large number of observations still has statistical 
significance. 

Obviously both the uncertainty in modelling the nongravi- 
tational forces as well as the bad determination and separation 
of the parameters, which presently still are the main source of 
uncertainty in the prediction of the perihel time for 1986, 
will not cause much lore essential positional uncertainty over 
shorter time intervals in 1985-6 (see also chapter 2.1). If 
light biases definitely do not exist then one could expect a 
position accuracy of 100 ... 150 km. 

The main source of position uncertainty, however, will be 
the possibility of a light shift of unknown magnitude, as 
suggested by the results of cases 4 to 6. If it cannot be taken 
into account explicitely, one must expect uncertainties up to 
700 km. As is visible from comparison of cases 5 and 6, the 
accuracy can improve essentially by increased observational 
effort from the end of 1985 to the first part of March 1986, 
down to below 200 km uncertainty. 

Even if the exclusive use of the position angles appears 
to be the only possibility to realistically guess the amount 
of uncertainty in the case of a light shift of unknown lapse 
for error estimations, it remains to question whether this 
method should indeed be used on orbit computations in 1986. 
This question cannot be answered presently because it strongly 
depends on our knowledge of the essence of the light offset 
and the observations obtained until March 1986. If the amount 
of the offset is ascertained to be smaller, or is modelled 
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better than the calculated position uncertainty using all 
observational information (including the elongations), then 
this is good so. But if this is not the case the suggested 
method should be used, because the primary intention is not 
to decrease the formal position uncertainty but is rather to 
exclude systematic errors larger or even of order than it. 
Also, if the light offset does not coincide accurately with 
the radial direction then probably at least its main part 
and the corresponding part which causes the most uncertainty 
in the target plane will be eliminated, whilst tangential parts 
mainly influence the arrival time. Most desirable, however, 
is to obtain models for the light shift from the theory of 
cometary comae and thus a decrease of uncertainty similar 
to the cases 1-3 discussed above. 

The presented assignation of the light shift as the main 
source of position uncertainties is in certain contradiction 
to the results of D.K.Yeomans et ai. (1982). These authors 
adopted ad hoc a certain magnitude S for the light shift 
depending only on the heliocentric distance (200 and 1000 km 
at r = 1 A.U., varying to r" 2 and r~ 3 , respective), 
subtracted it from the simulated observations, fitted through 
an orbit without solving for S„, and represented the result 
for the corresponding position uncertainty for March 1986. 
Subsequently it was concluded, that the light shift has 
little influence on the position accuracy in March 1986, because 
only the projection of the true uncertainty in the line of sight 
(which then coincide approximately with the radius vector) 
enters into the observed position. However, such investigations 
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only adress the problem in a very limited way because the 
assumption of, and the fitting of the observations with, 
a certain lapse of the light shift causes essential 
systematic errors not canceled on the fit and not indicated 
by the somewhat increased mean errors (in a simulation under 
similar conditions appearing to approximate the half magnitude 
of S„). Thus, the effect of the light shift to the geocentric 
position in March 1986 has only little to do with that to the 
heliocentric position. 

The error ellipsoides given by these authors for the cases 
without light shift agree in whole with ours in orientation 
and the ratios of the semi axis, but are a few times larger. 
This is probably caused by more pessimistic expections for 
the observations. 

Using observations 1835 to February 1984, r.A.Morley (1984) 
gave for the corresponding determination of the position in 
March 1986, A = 5700 km, B = 54 km and C = 15100 km, 
which is not in good agreement but of the same magnitude as 
our results to case Ob in table 12. 

Additionally, it was interesting to check how far the 
position accuracy can increase by special ob< 



from end of 1985 until March 1986. 



jservational effort 
One possibility for 
such, according to a suggestion by E.Boweii at Lowell 
could be the observation of 



Observatory, 
appearent close encounters of comet 
Halley to background stars. If these stars previously have 
been observed by transit circles, and if star-like images of 
the comet have been obtained by very short exoosures, a position 
accuracy of 0V1 is possible. For the error estimations, 
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from the middle of October until the end of 1935, 80 appearent 
encounters were assumed, for March 1986 further six encounters, 
of which each follows a position to ± V 1 (the number of stars 
appropriate by its position and magnitude difference to the 
comet is much larger). The other assumptions correspond to 
those of case 1 and 6 above; in particular, in case 6 for only 
the first twenty encounters the elongations to sun were also 
used. The results (no. 7 in table 12) correspond approximately 
to the expections regarding the increment of observation weights. 
We do not intend to enter into the technical particulars of such 
an observation project and some limitations in obtainable 
accuracy, although consideration of these would not change 
the result of a considerable improvement of the position accuracy 
by such a project. 

A further question was, whether it is imperative to abandon 
completely the observations of earlier apparitions and use only 
those of the present one. 

The results in table 13 are given only for the best and worst 
assumed cases. In the latter one, two assumptions about the 
observations in March 1986 were made to see their influence. 
Now, some of the parameters should not be taken as unknowns 
but should be considered with certain mean errors instead. 

The comparison between the results of the four first cases 
in table 3, in which always the six orbital elements had to be 
estimated, again suggests strongly that in 1986 the nongravi- 
tational forces will have only very minute influence on the 
position uncertainty. Also an uncertainty of any perpendicular 
force parameter A 3 as large as four times its present uncertainty 
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will cause only an increase of the uncertainty in z-direction 
by 15 km in the favourable observational case, but is completely 
succombed by other uncertainties in the less favourable cases. 

By comparison of these results with those in table 12 
it is evident that, by not taking into account the earlier 
apparitions, the uncertainty is increased two to five times. 
If values for both the nongravi tational parameters and the 
revolution time are considered instead of estimated, the 
increase is only 1.5 times. However, even if one solves 
only for one unknown, the perihelion time, under consideration 
of all other parameters as known without any uncertainty, 
the position uncertainty is not significantly below the limits 
attainable by a general solution including the earlier 
apparitions and increased observational efforts. Because, 
furthermore, by considering values and mean errors of several 
unknowns one would fall back upon the earlier observations 
implicitely, and because without considering a part of the 
unknowns, the computed uncertainty is much larger than the 
systematical errors produced by fitting observations from 
different apparitions using the available force models, 
it appears not to be imperative to use observations of the 
present apparition exclusively, but it could become appropriate 
later not to use the apparitions before 1835 or 1759 (cf. 
chapter 2.1}. 

5.2. Systematic errors in the predicted positions 

Although systematic errors are not the object of an error 
computation, we shall at this point make some remarks as, in 
few cases, one can guess their amount. 
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5.2.1. Different models for the nongravit at ional forces 

To get an idea about the influence on the position of the 
comet due the errors of assumed force models, one could 
compare the results after application of different models 
to one and the same observations. In particular, these 
differences depend completely on model and observation 
distribution, much more than the formal position uncertainty 
considered above. Neverthless, some results will be given 
as examples to show the order of magnitude of these differences. 

The differences in position predictions which resulted from 
a fit using two theories according to H.Rickman and c.Froeschle 
(1982) are given in table 14 and are compared with the results 
using the sublimation formula by Deisemme which was also used 
for the simulation of the expected observations. The assumed 
observational distribution approximately corresponds to no.l 
on table 11, but in the second case of table 14, after October 
1985 only position angles were used. These results suggest 
that the lack of knowledge of the accurate force lapse in 1986 
will produce only minute position uncertainties. Of special 
interest is the strong dependence of the calculated position 
uncertainty on the adopted model. This is mainly due to the 
different lapse of the parameters in the three models, wherein 
they are determined by very different portion of the observations 
(e.g., for Rickman-Froeschle with thermal inertia of the nucleus 
1^^=130, much more by the far distant observations then on the 
other models), and also by very different correlation with the 
other unknowns. 
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Also, using a model which is sufficiently accurate 
over the heliocentric distances covered by observations 
and applying this to the apparitions 1607 - 1986 and 1835 
- 1986, respectively, will give some differences in the two 
results for perihelion time in 1986 and the other parameters. 
These differences depend mainly on the error of assumptions 
about the secular behaviour of the parameters, e.g. in the 
above cases we assumed linear dependence (1607 - 1986) or 
parameters constant with time (1835 - 1986), respectively. 
It is not presently possible to conclude anything about the 
amount of the corresponding position difference in 1986, which 
is essentially A_r = jr-iT + i£ (r position of the comet in 
March 1986, aT,a_tt differences of both results for the peri- 
helion time and location of the perihel). The corresponding 
differences for 1835 and 1910 (using observations until 1984) 
let us conclude that this can amount to a few hundreds of km, 
but this we must see in 1986. Whilst presently solutions 
1607 - 1984 for the predictions appear preferable, it could 
become possible to prefer solutions 1682 or 1835 - 1986 then 
because of the better elimination of time dependence of the 
force parameters and the better fit of the observations 
although these would only yield a minute increase in the 
calculated uncertainty. 

Because of short-term fluctuations of the nongravitational 
forces we do not expect position shifts to be important. 
The magnitude of the nongravitational forces is approximately 
]A| = 0.1 km/d 2 . Even if a perturbation of the same magnitude 
as the main force is acting always in the same direction, after 
one month it would have produced a position error of only 
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-v-100 km. Such long-term effects, however, would have to 
be considered as a part of the (global) force theory, and 
furthermore, by the continuous fit to the observations 
they would produce only much smaller errors in the position 
predictions, as discussed above (cf. table 14). 

D.K.Yeomans et ai . (1982) have investigated extensively 
the positional effects of fluctuations composed by a decaying 
earlier fluctuation and a new random one. An amplitude of 
20% of the main force and a time scale of one day was choosen. 
However, considering the above magnitude of the nongravi tational 
forces on comet Halley, one can already compute by head that 
such perturbations are only of subki 1 ometri c amounts, so that 
similar computations were not made. 

5.2.2. Different reference systems for observations and 
coordinates of earth 

A source of essential systematic errors in the use of 
different reference systems for observations and earth 
coordinates. This is of special importance, because the 
navigation of space probes is bound to the earth rotation, 
so that the position of the comet must be known with reference 
to the dynamic equator and equinox and errors do not cancel 
but enter on the targeting accuracy. To the present state, 
the FK 5 system coincides with the dynamic reference system 
better than to O'.'Ol and thus is sufficient for our purpose. 

If for both the observations and the earth coordinates 
the FK 4 system, for example, is used, the error in the 
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position values corresponds to the transformation FKk-FKs, 
although small dynamical inconsistency exists because of 
the inertial rotation of the FK^ system. In practice, 
however, bases refered to different systems, e.g. observations 
in the rotating FK^ system and the approximate inertial resting 
earth coordinates according to new radar theories (0E102, 
DEI 18 etc) are used. Then the corresponding errors in the 
cometary position are no longer independent from the 
observations and can be estimated only by simulation 
computations . 

Subsequently the case was considered where the positions 
of earth corresponds to the theory DE118, whilst the observations 
of the comet refers to the FK^ system but are corrected for 
the elliptic aberration. It was noted that such bases are 
often in use in practice. The calculated positions for the 
comet were compared with corresponding results, assuming that 
the earth coordinates are refered to the FK5 system, and using 
the same observations of the comet reduced to the FK5 system, 
too (by application of the correction of equinox, of precession, 
of elliptic aberration, and of approximate local corrections). 

The results for four different assumptions of the observation 
distribution are presented in table 15. For essential 
simplification, the differences explained above were not 
taken into account on the observations of earlier apparitions 
again. Certainly, they do not enter into the result for the 
perihelion time for 1986 but do enter into the other elements, 
especially into the orbital plane orientation, which is of 
relevance in this case. Insofar, these results may have an 
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error of 20% or so. However, they are sufficient to show 
that here one has to expect considerable systematic errors 
whose accurate values cannot be taken into account other than 
by direct consideration. Thus, it is urgently necessary 
to eliminate these discrepancies in the bases. 

This point is also of relevance for the later investigation 
of the nongravitational force lapse by means of positional 
observations. A correction by some V 1 in far distant 
observations corresponds to an considerable correction of the 
mean anomaly. 
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Date 



1982 10 


18 


4 


4 




4 






11 


16 


2 


2 




2 






12 


11 


2 


2 




2 






1983 1 


14 


1 


I 




1 






2 


13 


2 


2 




2 






1984 1 


29 


20 


20 




20 






2 


28 


1 


1 




1 






10 


29 


10 












11 


25 


20 












12 


22 


60 


10 




10 






1985 1 


18 


40 


20 




20 






2 


14 


40 


20 




20 






3 


13 


20 


10 




10 






4 


10 


4 












8 


4 


10 












8 


25 


60 


20 


20 


20 


20 


20 


9 


22 


200 


40 


40 


40 


40 


40 


10 


19 


200 


40 


40 


40 


40 


40 


11 


14 


200 


40 


40 


(40) 


(40) 


(40) 


12 


7 


200 


40 


40 


(40) 


(40) 


(40) 


1986 1 


1 


80 


20 


20 


(20) 


(20) 


(20) 


1 


10 


20 












3 


5 


10 


2 


20 


(2) 


(10) 


(2) 


3 


8 


10 


2 


2 


(2) 


(10) 


(2) 


3 


11 


10 


2 


2 


(2) 


(10) 


(2) 



Table 11 — Assumed observation distribution until March 1986 



Error estimations of the position on 1986 March 13 are provided 
for six different assumptions about observation distribution until 
then. The table gives the assumed distributions of observation 
weights. One weight unit corresponds to one observation with a 
mean error of l'.'O . 
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Taile 12 — Mean error (la-variance) of the position predictions 
for comet Halley in March 1986 



Obs . : assumed observation distribution until March 1986 (see table 11) 
a: included observations 1607 - 1911 and solved for 

9 unknown parameters (elements, Ai,A2,S) 
b: included observations 1835 - 1911; 8 unknowns only (without B) 
c: included observations 1607 - 1911 and 80 assumed observations 
for mid October until end December 1985 and 6 further for 
first part of March 1986 with a mean error of O'.'l each 
Case no. refers to the observations which have been presented until 
February 1984. For comparison, on top are given the results 
from the observations 1607 to 1911 only. 
aX,ay ,CTZ,ar: mean errors of the heliocentric equatoreal coordinates 

and distance of the comet at 1986 March 13.60 UT 
A,B,0: semimajor and semiminor axis of the error ellipse in the target 
plane of 'Giotto', and direction of the semimajor axis 
(ecliptic 0° , orbital plane of comet Halley 4") 
C: mean positional error of the comet in flight direction of 'Giotto' 
oT,aq,ae,ou,an,ai : mean errors of the orbital elements of comet Halley 

at osculation epoch 1986 February 19 
oa,afi,aA: mean errors of the equatoreal geocentric coordinates of 
comet Halley at 1986 March 13.60 UT 

See also table 16. 



Table 13 — Mean error (la-variance) of the position predictions 
for comet Halley without use of the observations at 
earlier apparitions 

Corresponding to table 2. In the last column are given the parameters 
which have been considered with certain mean errors, instead of estimated. 
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No. AX Ay Az aX ay aZ 



1 











276 


591 


96 


2 


-3 


+24 


+6 


274 


549 


85 


3 


+2 


+6 


+2 


283 


587 


93 


r 











1737 


2852 


1168 


V 


-11 


-17 


-7 


840 


1140 


462 


V 


+20 


+20 


+ 11 


1423 


2174 


905 



Table 14 — Systematic differences in position due to application of 
different models for the nongravitational forces 

For the models by Delsemme and Rickman-Fzoeschle (I =130 and 1000) 
are given fNo. = 1 ,2,3, resp.): 

AX.Ay.Az: positional difference to the results using the model by Delsemme 
oX,ay s az: mean error of the position 

The first three line refers to computations without, the second ones 
with considering a possibly present light shift between nucleus and light 
center of the comet, by using for observations at r<2.0 AU only the position 
angles to sun in the latter case. Moreover the observation assumptions 
corresponds approximately to case no.l in table 1. 
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No. AX Ay AZ Aa A6 (as | |At] 

km km km km km 

1 +370 -240 +30 -O'.M -O'.'l 410 170 

3 +280 -150 +23 -0.0 -0.0 290 140 

5 +560 -160 +90 -0.1 +0.4 470 350 

6 +1040 -204 +190 -0.3 +1.0 810 710 

Earth -48 -380 



Table 15 — Systematic errors in the position predictions due to using 
observations refered to the FK4 system and earth 
coordinates according to DEI 18 

Presented are the differences in sense a)-b) of position predictions 
based on; a) coordinates of earth refered to the FK 5 system, cometary 
observations corrected for the correction of equinox (FK 5 -FKi+), correction 
of precession (IAU 1979 - Net/comb), and for the approximate local 
corrections FK 5 -FK H , b) earth coordinates corresponding to DE118, 
observations not corrected. 

AX.Ay.AZ: differences in the heliocentric equatoreal coordinates of 

comet Halley at 1986 March 12.60 UT, refered to mean equinox 
1950 

Aa,A6: differences in the residuals of observations at 1986 March 13 

AS.At: difference in flight direction and target plane of 'Giotto' 

For comparison, the position difference FK 5 -DE118 of earth is given 
at the bottom. 
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a) orbital elements (epoch 1986 Feb 19.0) 



No. T(TDB) 



1 1936 Feb 9.50954 0.5871050 0.96727924 111°84726 58°. 14370 162°.23921 

2 1986 Feb 9.46133 0.5871017 0.96727381 111.84662 58.14388 162.23943 



No. A A B 

1 -0.0172 +0.0084 0.015973 ±0.000014 -0.00502 ±0.00016 

2 +0.1110 ±0.0192 0.015564 ±0.000021 



observations 
No. arc no. a y 



1 1607 - 1984 111 1'.'49 1!!25 S FK5 ITA 

2 1835 - 1984 104 1.01 0.88 S FK5 ITA 



b) mean errors (referred to ]i-l "O) 



aT 


6148-10" 6 d 


aii) 


1594-10 -7 " 


OX 


22390 km 


A 


17200 km 


aq 


771-10-9 


a£2 


897-10" 7 ° 


oy 


2180 km 


B 


60 km 


ae 


91-10" 3 


ai 


481-10 -7 ° 


oz 


5650 km 


C 


15540 km 


aAi 


0.0218 






or 


13520 km 


e 


20° 


aA2 


2.349-10" 5 















Table 16 — Updated orbital elements and error estimation 

The orbital elements correspond to those of no. 1 and 3 in table 3, 
the mean errors to the 1835 - 1984 solution and no. Ob in table 12. 
In comparison to the earlier solutions, few observations obtained 
recently have been added, but the weights of some observations made 
in early 1984 have been decreased, and the declinations of the 
remeasured 1911 observations (cf. chapter 1) were omitted. 
Due to this, the weight sum of the observations decreased by 7%, 
and the actual mean errors (tabulated values times ]i) have not 
become smaller. 
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Appendix. A program for the computation of the mean error 
(la - variance) of position estimations for comet Halley 
at the Giotto encounter 



I. Introductionary remarks 

The calculation of the mean error of parameters determined by application 
of the Method of Least Squares, and of functions of these unknowns, can, 
for example, be obtained in the following basical wav (e.g. for technical 
perforation, see G.ueutier, 1982). From the system of normal equations 
for the unknown parameters u_, 

1 = S"ii a) 

the variance-covarianee matrix ^, correlation matrix K, and the mean errors 
c of the m parameters follows by 

V.. = ai K..a j = ^g- I ) ij b) 

with K. -=1. The mean error of an unit weight, u, can be guessed by 
2 s a, s° - 1-u. x 

U = TS =^= C) 

n-m n-m ; 

(n number of equations of condition, s°,s residual squares sum before and 
after the orbit improvement, respectively). By substitution of one of the 
unknowns by a function f = f(u) of them, the transformation of C, V and a 
results in 
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m m , - 
ff f ^^P^^iVjsu. d ) 

which permits the computation of the mean error of a function of 
the unknowns. 

For the calculation of the mean error of ephemeris places obtained 
after an orbit improvement, the partial s of the position values to the 
unknowns and the normal equations of the orbit improvement are need. 
Thus, the only suitable place to apply a corresponding programme is 
to include it to the orbit determination program and call it subsequent 
to the orbit improvement(s). For this reason this is assumed to be the 
case on the programme explained below. Moreover, it was assumed that 
for usual orbit determination problems (e.g., the computation of the mean 
errors of the elements if using rectangular initial values as parameters) 
the user already has a subprogramme for computing the mean error of a 
function of the parameters, or otherwhise can write this quickly by using 
formula d) above. This subprogram subsequently is refered as MF2. Then 
the enclosed programme, called ERAN, contains only the calculations 
necessary in specific regard to the position uncertainty of comet Hal ley 
on 1986 March 13.6 UT, the assumed time of the Giotto encounter. 

II. Explanation of program ERAN 

The mean errors of the heliocentric equatoreal coordinates x of the 
comet, immediately follow by application of MF2 to 3x/au_, the partial s 
between position and parameters u. The latter have to be computed by 
the orbit programme in the same way as the equations of condition for an 
observation at this time, and transfer to ERAN; however, they are already 
predetermined ^jery precisely and could instead be taken from the example 
given below (assuming elliptic elements as the parameters). 
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For computation of the mean error of the heliocentric distance r, 
the 3r/3£ follows immediately by 

3£ * 3£* 3X_ ' 

where sr/ax^ = (x,y,z)/r for the cometary position is also pre-determined. 

Let be x':= (r,s,t) the position of the comet in a coordinate system 
in which the flight direction of Giotto is perpendicular to the r/s-plane 
(target plane) and parallel to the t-axis. Then the mean errors of x_' 
can be computed by application of MF2 to 

M'-M.M' -n 

3£ ~ 3£3X ' 

The orientation of r,s (miss vector) in the target plane is arbitrarly. 
In prcatice, the position uncertainty in the target plane has to be 
computed for different directions of the miss vector. 

If the flight direction of Giotto is called y_ = (v cos a cos 5, 
v sin a cos s , v sin &), then the transformation x -*■ x_' is obtainable 
by rotating first the z-axis by +a, and then the new y-axis by 90° -6. 
After this, the new z-axis coincides with the flight direction, and all 
orientations of the miss vector in the target plane can be obtained due 
to additional rotation of it by an arbitrary angle e running from 0° to 
360". Thus, we have 



/ cos sin o 0\ / sin & -cos fiWcos a sin a 0> 

4^ = f-sin o cos e 1 -sin a cos a | g) 

\ 1 ) \ cos 6 sin s/\ 1 



3X 



?• I 



The uncertainty in t-direction is independent on e, and if e runs 

from 0° to 180°, only the mean error in the r-direction must be computed. 

The angle a(s) between the corresponding r-direction e and a certain 
equatoreal direction .n (ij]J = 1) can be computed by 

cos a(9) = el-n 1 =(0| • |f'n = rC|| { n. h) 
-r- ^oj 3 - i=l 3x i 1 

If £ is the normal vector to the ecliptical plane, then at cos a(e) - 

the r-direction lies in this plane, and the corresponding r-uncertainty is 

that of the miss-vector along the direction of the node between ecliptic 

and the target plane. 

These are essentially the formulas after which ERAN works. A listing 

of the program is enclosed, as well as a table which give the used notation 

and references of the variables used in the program. 

III. Example 

As an example which can be used to check out ERAN before application 
(and to check a procedure for formulas a),b},d), if necessary), are enclosed 
the values for the normal equations, correlation coefficients, mean errors of 
the parameters, and results printed by ERAN, corresponding to case no. lb 
on table 12. Note that all results are normed by setting y=l "0=1/206264. 8, 
instead application of formula c) . The parameters are the elliptic orbital 
elements and two nongravitational parameters, expressed in the common units. 
For example, the correlation coefficient between e and k x is -0.9977339, and 
the partials of the cometary position to the perihelion time corresponds to 
the negative velocity of the comet. Although not necessary, ERAN prints the 
uncertainty in miss vector until 0=360° . The node between target plane and 
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ecliptic (used as reference point for the orientation of the error 
ellipses on table 12 and 13) is at 6=107°, the semimajor axis at 6=162°. 
The semimajor and semiminor axes amount to 54 km and 21 km, respectively, 
and the uncertainty in flight direction 123 km. 
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Variables in subprogram ERAIJ (cont.): reference listing 
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Listing of subprogram 



Subroutine fran 

1966 SUBROUTINE ERAN KALLEY 

1967 IMPLICIT DOUBLE PRECISION (A-H,0-Z> 

1968 C Berechnet die Un s i c her hei t en der Pos it ions uer t e des. Hi m*t I skturwrs. 

1969 C Drojeziert in eine vorgegebene Richtunq CKomnonent en y im Aeouator- 

1970 C system ausgedrueck t ># so dass diese Rictituna mr t-Achse eines neuen Koordi- 

1971 C m tensyst ernes (R»r#s>t) und die Ausriehtung dessen r-Achse variiert uird 

1972 C (Drehuinkel W* entsor, Neigung der r-Achse 2U einer vorgeaenenen Ebene 
■J 973 C uie etua Ekliptik oder Kometenbahn, 03LI1). 

1974 C Seispiel: Beim Vorbeiflui einer Sonde an einem Xoi»eten spLl die L'nsi- 

1975 C cherheit seiner Position in F luar ichtuna ( t-R i chtung) und in 

1976 C miss vector <bei ?aehlunq von Ekliotik Oder Koine t enbahnebene > anje- 

1977 C qeben verrten. 

1973 C SD,CD,SR,CR Sinus und Kosinus der Rfkt. und D e VI . in ft ichtuna Xr 
1979 C X1» X2 Norma Ifnvktor senkrecht auf Ekliptik bzw. Komet enbahnefcene 
1R8Q C i m Aeauat orsys t em Cueitere Ebenen entsprechend einbaubar) 

1981 PARA»ETFR (LU*1 1 > 

1952 C01M0N/ERRAN/ , WEIT/DXDtjaLI»3> 

1953 DIMENSION X(3),e{3)* X1(3>»X2(3).-FR (3) , FX (31/ DRDX(3,3) ,D° O'.J CLU/3) 
1984 DIMENSION X3(3) ,DDD|I(IU) ,Dt)0X(3> 

19S5 DATA Z'JEIPI.AE 16. 2 831 3530 71 79586C3 » 1 45 597s 70.65255 on/ 

1986 DATA X /♦.69821F071978D0,-. 66212925993100, -.272lfc9715Q64D0/ 

1987 DATA XI /.0D0,-.3978 l i1153349Dr,+.91743S9645(11D0/ 

1988 DATA X?. / + .Z5?3?4')09205nP»+.2312127ri537SDn»-.93776«!9016«7DG/ 

1989 DATA X3 /-.028760S32565DR.-.40579?4S57S3DO,+ .?1S509491953I)C/ 

1990 DATA DDDX /- . 5 8S3434 DOr- . 73 221 3730, - .343 0?S P. 00/ 

1991 C (Entsoreehen der Vorbe i f I uqr i c h tuna »cn Giotto an Haltey'schen Kerne ten 

1992 C und der Normalen senkrecht zur Ekliptik, HaUey- und S i ott obahne bene) . 

1993 C Bezoaen auf mittC. Fehter von 1"s 

1994 AE=AE/206264.8DQ 

1995 IFOI.EQ.rnBETUSN 

1996 SP=X(3> 

1997 CO=DSQRT(1.D0-SD*SD> 

1998 SR=X(2)/CD 

1999 C CR=X(1>/CD 

2000 CR=0SORT(1 .D0-SR*SR) 

2001 IF(X(1).LT..0D3)CR=-CR 
22C? URITE(ft,12)ZEIT,M 

2003 12 FORMATC/' Ab I ei tungs mat r i x des Ones zur Zeit ',F13.4» 

2004 . ' nich den ',12,'' Unbekann t en : ■ ) 

2005 D013L=1»3 

2006 13 URITE(1.U>(r,XiiJ<I#Ll»I = 1»f ) 

2007 14 F0PMATt7C1X.G17.10)) 

2008 C "ittl. c ? >,ier FX in aequat orealen Koordinaten X,Y,Z»» des Kometen 

2009 DO10J=1,K 

2010 10 B0DU(J)=nXDIIU,1 )*Dr>DX(1 )+0X0'J(J,2) «DD!>X (2>*D¥DU<J»3> *00DX C35 

2011 CALL »F2(",DDDy,FD) 

2012 6021=1,3 

2013 2 CALL iF?<«!»»Xaun.I)»FXCI)> 

2014 URITE(6,20)ZEn,FXC1>*AE,FX(?>*AE,FX(3>*AE»FD*AE 

2015 20 FORMAT'///' 'Jns i cher neit der aeaua t ore al en Xoordiniten ur1 ', 

2016 . 'hetioz.Distanz: ',11' Ze i t: ' ,F1 5. 3, ' <i* ',FS.1,' dy '» 

2017 . FS.1.' dz ',F8.1,' dr ',F8.1." (in km)'/) 
2013 URITE(6,21) 

2019 21 F0Rf«AT(/72t^Xr'F•hlereUip5oid , //' Wi nk el ' ,5X» 'dr ( km) • ,5X , 

2020 . ' dt (km)" ,6X»'C1 , .7X. , C2' ,7X,*C3' /) 

2021 C Serechnuno. von DP3X (I ,J)=dR( J >/dX (I) "nit Variation von W 

2022 C (ergibt die Fehl ere 11 i ose in der s»t-Ehen-.) 

2023 DRDXd,3) = CD*CR 

2024 0R0X(2»3)'C5*5R 

2025 DRDX(3»3)=SD 

2026 DOHU»1»36 

?027 W»ZWEIPI/36.D0*IW 

2028 SJ^DSIN(U) 

2029 CW=DC05(u) 



79 



Listing of subprogram ERAN (cont.) 



2030 DRDXC1/1)=CU*SD*CR-SV*SR 

2031 DRDXC2.1 >»CW*St>*SR+SW*CR 

2032 DRDX«3..1)=-CW*CD 

2033 C 0R&XC1.2>=-SU*SO»CR-Cy*5R 

2034 C DRDXC?,2)=-SV*5D*SR+C«»CR 
?ni5 C DRDX(3»2)=SW*Cn 

2036 C Konrrolle: bei Transformation der Richtunq X» wird t = 1*s=r-0 

2C.37 C 1)091 = 1,3 

2035 C R(n = .TDO 

2039 C D0'J = 1.,3 

2040 C 9 R(I 1=RC mX(J)*DRt)X< J*I) 

2041 C = erechnune> des iiinkels 9I.tr der r-Achs* ?u den voraet). Ebrrrn (Ekliotik* 
2^42 C "Comet eibahn usu). Falls cos <0BLO ) =r- Hoipoonent e d. Notts I en=C lieat 

2343 C di e r-Achse in der b#tr. Enene 

?C« S0DLQ1 = >'1<l>*nR7X(1,1)+X1C2)»ClRDXC2, T) +X 1 <"=:) *DRDXC3* 1 3 

2 04 5 S0RLO2=i'2(1)*DS!>y(1»13 + y2(2)*DRDX(2,1)+X2r3)«DRDX{3,15 

2046 SC9LQ , ! = X3C1)*0PnX(1/1) + X3(?)*PRDX{2/1)+X3C3)*l)ft0X(3,1) 

2047 C Tr »nsfortnti onsmat r i« PDU ( I> J ) = dR C J ) / dM (1 ) , fyer r immer* fuer s nicht/ 
2^4? C fuer t nur einmal ( be i IW=1) zu berecinen 

204 9 J=-l 

2050 7 J = J + 2 

2051 D06I = 1,« 

2D52 DR5UCI, J) = .Ql>rt 

2053 DP6k>1e3 

2^54 6 DP5UtI»J> = t>RI>UU,J)+DR3X(K.J>*5XDUCI /if) 

2055 C Fehter in r-Richtunq 

2055 CALL «F2CM#DRtHJ<1,J>,FRU)> 

205 7 IFCIW.EQ.1 .A-JD.J.EQ.I 1G0T07 
2053 C (usdruek der Erqebnisse 

2 05 9 WRXTEC*.22)10.*IU»FRl1)*AE,FR(3)**E#S0BL«i.SSBLQ2..SCBLQ3 

2050 22 F0R,»AT<2X,F5.1/2<3X/FS.1),3C2X/F7.3)) 

2061 1 CONTINUE 

2052 WRITEC4/23) 

2063 23 F0R»»T<//' dr,dt mittl. Fehler in r- b?g. t-Ei c htung, ■ , 

2064 . /' C1-C3 Sinus der Neiqung der r- Fi chturtq zu Ekliptik'» 

2065 . ' bj w . Kometen- und Gi o ttohahn ' // ) 
2065 RFTURN 

2057 C DE3US SU?CHK 
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